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FOREWORD 

This  report  describes  an  in-house  effort  conducted  by  personnel  of 
the  Microwave  Technology  Branch  (AADM),  Electronic  Technology  Division 
(AAD),  Avionics  Laboratory  (AA),  Air  Force  Wright  Aeronautical 
Laboratories,  Wright-Patterson  Air  Force  Base,  Ohio,  under  Project  2002, 
"Microwave  Technology,"  Task  200202 ,  "Thermionics,"  Work  Unit  200202 44, 
"Physics  of  Matrix  Cathodes." 

The  work  reported  herein  was  performed  during  the  period  1  October 
1976  to  30  December  1980,  under  the  direction  of  the  author,  Capt  Jess  B. 
Scott  (AFWAL/AADM-1 ) ,  project  engineer.  The  report  was  submitted  by 
the  author  in  January  1981. 

The  author  is  indebted  to  Dr.  G.  K.  Medicus  for  helpful  discussions 
and  guidance  during  the  course  of  the  theoretical  work,  to  Lt  D.  L.  Allen 
for  assistance  in  maintaining  the  computing  equipment,  and  to  Lt  J.  Targove 
for  assistance  in  some  of  the  programming. 
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SECTION  I 
INTRODUCTION 

One  purpose  of  thermionic-emission  theory  and  space-charge  theory 
is  to  relate  in  a  sensible  way  the  current  drawn  from  a  cathode  to  the 
voltage  applied  to  the  anode.  Tests  are  made  on  cathodes  in  geometries 
kept  as  simple  as  possible  to  make  it  easier  to  relate  measured  current 
density  J  as  a  function  of  voltage  V  back  to  the  basic  characteristics  of 
the  cathode.  The  usual  cathode  property  of  interest  is  the  cathode  work 
function  or  the  zero-field  current  density  Jq  which  is  an  exponential 
function  of  *  and  the  temperature.  Once  the  basic  properties  of  the 
cathode  are  known  the  performance  of  the  cathode  can  be  predicted  for 
the  desired  application. 

Different  techniques  are  available  for  relating  the  zero-field 
current  density  JQ  of  the  cathode  to  the  J{ V }  measurements  depending  on  the 
magnitude  of  J  and  the  range  of  voltage  used  in  the  test.  The  develop¬ 
ment  of  theory  covering  three  modes  of  operation  -  two  which  neglect  the 
effect  of  space  charge  and  one  which  includes  space-charge  effects  -  is 
reviewed  below.  (1)  If  the  voltage  is  sufficiently  large  to  overcome 
space-charge  effects  in  the  interelectrode  space,  the  current  becomes  an 
exponential  function  of  the  square  root  of  the  voltage.  This  relationship 
between  the  current  and  voltage  is  known  as  the  Schottky  effect.  (2)  If 
JQ  is  sufficiently  small,  space  charge  effects  can  be  neglected  for 
values  of  J  up  to  JQ.  For  this  experimental  condition  the  data  do  not 
show  the  increase  of  current  due  to  the  Schottky  effect  until  the  voltage 
is  increased  far  beyond  the  value  needed  to  obtain  JQ.  In  this  case  the 
current  is  an  exponential  function  of  the  voltage  up  to  and  then  a 
very  sharp  saturation  for  J  greater  than  JQ.  (3)  If  space-charge  effects 
are  significant  at  JQ,  the  current  departs  from  the  exponential  behavior 
at  low  voltages  and  rises  to  J0  according  to  the  space-charge  theory. 
Following  a  review  of  the  historical  development  of  theory  for  these 
three  modes  of  operation  is  an  overview  of  the  new  additions  to  this  theory. 
The  new  additions  to  the  theory  are  then  developed  in  detail  in  the  body 
of  the  report. 
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1 .  BACKGROUND 

The  exponential  nature  of  the  current  as  a  function  of  voltage  in  the 
retarding  field  range  is  discussed  by  Davisson  (Reference  1).  If  space- 
charge  effects  are  small,  the  current  density  saturates  abruptly  and 
remains  steady  at  J  =  J  .  Shelton  (Reference  2)  gives  experimental  data 
which  fit  this  theory  very  precisely  and  shows  how  to  evaluate  the  basic 
cathode  properties  in  terms  of  the  experimental  data  for  his  particular 
example.  Medicus  (Reference  3)  generalizes  the  work  of  Shelton  and 
Davisson  to  develop  a  systematic  method  for  analyzing  retarding  field  data. 

Schottky  (Reference  4)  explains  the  exponential  rise  in  current  with 
the  square  root  of  the  voltage  for  large  voltages  (accelerating  field  at 
the  cathode).  In  the  absence  of  space-charge  effects  the  applied  potential 
interacts  with  the  surface  fields  near  the  cathode  to  lower  the  work 
function.  The  fact  that  the  surface  field  is  accurately  represented  by 
the  image  of  the  electron  in  the  conducting  surface  has  been  established 
by  the  agreement  of  Schottky  theory  with  experimental  data  (Reference  5, 
p  109).  Hutchison  (Reference  6,  p  305)  gives  a  review  of  Schottky  theory. 

Langmuir  (References  7  and  8)  summarizes  and  extends  previous  work  to 
give  a  theoretical  treatment  of  the  effects  of  space  charge  in  the 
interelectrode  space.  Langmuir  develops  functions  n  and  r  representing 
the  voltage  and  distance  in  the  space  between  a  paral lei -pi ane  cathode 
and  anode  and  outlines  a  series  of  steps  to  calculate  the  current  as  a 
function  of  voltage  for  this  type  of  diode.  Van  der  Ziel  (Reference  9) 
shows  that  Langmuir's  n,  >.  functions  can  be  extended  to  cover  retarding- 
field  and  accelerating-field  conditions  as  well  as  space-charge 
conditions,  however,  he  neglects  the  effect  of  the  electron  image  field. 

Van  der  Ziel  mentions  the  transition  point  from  retarding  field  to  space 
charge  conditions  and  names  it  the  "exponential  point".  Although  Van  der 
Ziel  gives  a  method  of  computing  n  and  '  for  accelerating  field  conditions, 
his  theory  still  predicts  that  the  current  density  is  identically  J  for 
accelerating  field  conditions  since  he  neglects  the  Schottky  effect. 

Kleynen  (Reference  10)  gives  several  useful  formulas  for  calculating  n 
and  ■'  and  extends  the  tables  of  r,,  •  values  given  by  Langmuir. 
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Ferris  (Reference  11)  takes  the  space-charge  theory  developed  to 

this  point  and  presents  it  as  a  family  of  normalized  JlV)  curves  which  can 

be  used  to  analyze  experimental  data.  The  experimental  data  are  presented 

in  normalized  form  and  superimposed  on  the  family  of  JM/1  curves  allowing  JQ 

to  be  determined.  The  normalizing  constant  introduced  by  Ferris  is  I,o 

which  represents  the  maximum  current  density  which  can  be  drawn  to  the 

anode  under  retarding-field  conditions  for  the  condition  that  I  tends  to 

0 

infinity.  Crowell  (Reference  12)  presents  a  slightly  amended  normal ization 
technique  which  simplifies  the  family  of  curves  presented  by  Ferris. 
Nottingham  (Reference  13,  p  65)  uses  I  ,  the  current  at  the  "exponential 
point",  as  a  normalizing  constant  to  give  a  family  of  JiV)  curves  only 
slightly  different  from  those  of  Ferris.  Nottingham  presents  experimental 
data  showing  that  the  second  derivative  of  J { V }  has  a  sharp  peak  at  the 
exponential  point.  However,  Nottingham  uses  the  second  derivative  only 
in  conjunction  with  his  family  of  JlV)  curves  to  evaluate  0Q.  Finally 
Rittner  (Reference  14)  presents  Langmuir’s  n,  f  tables  as  interpolation 
formulas  suitable  for  use  on  a  digital  computer  and  shows  that  Langmuir's 
theory  accurately  matches  experimental  J {VI  data. 

Aside  from  this  main  development  of  the  space-charge  theory  Crowell 

(Reference  15)  discusses  a  method  using  the  second  derivative  of  JIV1 

in  an  approximate  determination  of  J  .  Crowell  uses  a  numerical  study 

covering  a  wide  range  of  JQ  values  to  show  that  Langmuir's  theory  predicts 

an  inflection  point  in  the  J  •>  V  |  function.  This  gives  theoretical 

justification  to  the  standard  practice  of  associating  the  inflection  point 

of  experimental  J{VI  data  with  JQ  (IRE  standard  definition  of  inflection 

point  current  (Reference  16).  Crowell  gives  formulas  for  calculating 
2  2 

dJ/dV  and  d  J/dV  implying  a  functional  relationship  for  J(V}. 

However  Crowell  does  not  exhibit  J { V 1  in  functional  notation.  He  continues 
to  rely  on  a  series  of  descriptive  steps  for  calculating  J<V}. 

2.  OVERVIEW 

The  work  presented  in  this  report  extends  the  space-charge  theory 
of  paral lel-plane  diodes  in  several  ways.  Space-charge  effects  are 
included  in  the  accelerating  field  range  to  remove  the  incompatibility 
between  Langmuir  theory  and  Schottky  theory.  A  method  is  developed  to  use 
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a  measurement  of  Jg  at  the  exponential  point  to  calculate  without 
resorting  to  matching  a  family  of  curves.  A  number  of  properties  of 
Langmuir  space-charge  theory  is  described  in  much  greater  detail  than 
has  previously  been  done  through  the  use  of  functional  notation  to 
express  J •;  VI.  These  details  are  brought  out  in  a  discussion  of  the  fields 
and  charge  in  the  interelectrode  space.  As  a  side  light,  a  method  is 
developed  to  handle  the  edge  effects  in  a  class  of  paral lei -plane  diode 
testers . 

The  effects  of  space-charge  in  the  accelerating-field  range  are 
accounted  for  through  the  interaction  of  the  space-charge  field  with  the 
electron  image  field.  The  theory  calculates  a  J-V.  curve  which  rises  from 
J  =  with  the  Child's  law  slope  and  then  approaches  the  J • V 1  curve 
generated  from  Schottky  theory.  Since  Langmuir  space-charge  theory 
predicts  that  JiV,1  approaches  J  -  Jq  with  zero  slope,  this  theory  predicts 
a  discontinuity  in  the  slope  of  J  V  r  at  J  -  JQ.  This  discontinuity  is 
pronounced  in  the  J -  V  >  plots  for  low  values  of  Jo  while  it  occurs  over 
such  a  narrow  range  of  J  values  for  large  Jo  that  it  is  not  visible  in 
the  J;V-  plots.  However,  the  theoretical  discontinuity  shows  up  as  a 
sharp  minimum  in  experimental  second  derivative  plots  (Reference  13,  p  135) 
and  is  expected  to  be  apparent  even  for  large  JQ.  Observation  of  this 
discontinuity  should  aid  in  the  evaluation  of  cathodes  even  in  the  presence 
of  multiple  work- function  patches.  An  application  of  this  extended 
theory  is  given  to  show  how  a  common  "cathode  evaluation  technique"  lacks 
theoretical  justification.  The  theory  shows  how  this  "evaluation 
technique"  leads  to  erroneous  values  of  JQ  which  are  dependent  on  the 
diode  spacing.  (Reference  17  gives  experimental  data  demonstrating  this 
problem.)  A  close  examination  of  the  transition  point  from  the  retarding- 
field  range  into  the  accelerating  field  range  (the  "exponential  point") 

results  in  a  method  for  determining  J  from  a  measurement  of  the  current 

o 

density  at  this  transition  point.  This  evaluation  of  JQ  is  accomplished 
with  a  single  graph  of  JQ  as  a  function  of  Jg  without  the  matching  of 
experimental  J'V1  data  to  a  family  of  J T V 1  curves  as  was  previously 
necessary.  The  theory  calculates  a  discontinuity  in  the  second  derivative 
which  shows  up  as  a  sharp  peak  in  second  derivative  data  measured  by 
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Nottingham  (Reference  13,  p  135).  This  evaluation  technique  is  valid  for 
high-work-function  cathodes  and  can  be  applied  to  low-work-function 
cathodes  at  reduced  temperature. 

Previously,  Langmuir  space-charge  theory  has  been  used  to  calculate 
values  of  V  for  given  values  of  J  through  a  series  of  steps  (or  flow 
chart)  without  the  use  of  a  mathematical  expression  for  V  as  a  function 
of  J.  By  exhibiting  ViJ1  as  a  mathematical  expression,  trends  can  be 
observed  and  limits  and  derivatives  can  be  calculated.  As  an  example  of 
the  importance  of  functional  notation  the  braces  :  1  are  reserved  in  this 
report  for  designating  "function  of"  (for  example  ulO,  J  VI,  ln-xt). 

The  use  of  functional  notation  for  the  ViJ I  relationship  has  led  to 
the  description  of  a  number  of  properties  of  Langmuir  space-charge  theory 
in  much  greater  detail,  than  has  previously  been  done.  The  theory  shows 
mathematical ly  that  there  is  always  an  inflection  point  for  JiVi  in  the 

space-charge  range  and  this  inflection  point  coincides  with  a  discontinuity 

2  2 

in  the  second  derivative  d  J/dV  when  JQ  is  very  small.  The  theory  shows 
that  J  always  approaches  JQ  with  a  slope  dJ/dV  equal  to  zero.  The 
limiting  values  of  several  other  derivatives  of  current,  fields  and 
charge  density  functions  are  evaluated  at  the  exponential  point,  at  the 
transition  point  from  space-charge  to  accelerating  field  conditions,  and 
at  the  cathode  surface. 

As  a  side  light  the  report  gives  an  extension  of  Schottky  theory 
(without  space-charge)  to  include  the  effects  of  edges  in  a  class  of 
paral lei -plane  diodes.  The  field  solution  of  a  plane  edge  above  an 
infinite  plane  is  used  to  estimate  the  departure  from  paral lei -plane 
fields  near  the  cathode  edge.  For  a  typical  example  the  edge  effects 
introduce  a  66°C  error  in  temperature  calculations  based  on  Schottky 
theory  while  introducing  no  significant  error  in  the  determination  of  the 
zero  field  current  J  . 
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3.  THERMIONIC  EMISSION 


The  behavior  of  a  f ree-electron  gas  in  a  three-dimensional  potential 

well  is  sufficient  to  explain  the  basic  properties  of  thermionic  emission 

from  metallic  cathodes.  The  number  density  of  electrons  in  the  f ree- 

electron  gas  with  momentum  values  lying  in  the  interval  dp  dp  dp  is 

x  y  z 

given  by  Reference  6. 


N{p}  dp  dp  dp. 


...dPx  dVPA. 

exp  { [ (px2+py2+p2 


)/2m  -  Ep]/kT} 


:d 


where 

-31 

m  =  9.109  x  10  '  [kg]  =  mass  of  an  electron 
k  =  1.380  x  IQ'28  [J/°K]  =  Boltzmann's  constant 
h  =  6.625  x  10‘34  [J.sec]  =  Planck's  constant 


T  =  temperature  [ ’K] 

Ep  =  energy  of  the  Fermi  level  [J] 


Equation  1  is  the  product  of  the  number  of  states  existing  in  the  potential 
well  within  that  range  of  momentum,  times  the  probability  that  the  states 
are  occupied  at  temperature  T.  The  wave  nature  of  the  electron  is  used 
to  obtain  the  number  of  energy  states  available.  The  statistical  theory 
of  particles,  in  which  only  two  particles  are  allowed  in  an  energy  state 
at  one  time,  is  used  to  determine  which  of  the  possible  energy  states 
are  occupied  at  a  given,  nonzero  temperature.  The  Fermi  level  Ep  is  the 
highest  energy  state  which  is  filled  when  all  lower  energy  states  are 
filled  (at  0°K).  The  reference  level  for  zero  energy  is  the  bottom  of 
the  conduction  band  for  this  derivation. 


If  the  potential  well  has  a  finite  depth,  with  wall  height  <|>  electron 

volts  above  Ep,  Equation  1  shows  that  some  electrons  will  have  sufficient 

energy  to  escape.  For  these  escaping  electrons,  the  exponential  term  in 

-19 

Equation  1  will  be  exp  ie^/kTl  or  larger  where  e  (  =1.602  x  10  [C] )  is 

the  charge  of  an  electron.  For  thermionic  emitters  <}>  is  of  the  order  of  a 
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5 

few  electron  volts  so  the  exponential  term  is  of  the  order  of  1  x  10  . 
So,  for  the  escaping  electrons,  the  1  in  Equation  1  can  be  neglected 
compared  to  the  exponential  term.  Equation  1  then  becomes 

N{p}  dpx  dpy  dpz  -  (2/h3)exp{[EF  -  (px2+py2+pz2)/2m]/kT) dpxdpydpz  (2) 

For  convenience  let  the  normal  to  the  surface  lie  in  the  x  direction. 
Since  the  integral 

CD  O 

;  exp  {  -  Py  /2  mkT)  dpu  =  /2miikT  (3) 

for  u  =  x  or  y,  the  number  density  of  electrons  with  momentum  in  the 
interval  dpx  is  given  by 

N(pxldpx  =  (4-rmkT/h3)  exp{ Ep/kT}  exp  {-px2/2mkT}  dpx  (4) 

If  there  are  no  fields  external  to  the  potential  well,  electrons 

with  momentum  px  can  escape  where 
o 

p  /2m  =  Ep  +  e*  (5) 

o 

Then  the  total  current  escaping  from  the  potential  well  is  the  charge  of 

the  electron  e  times  the  number  density  of  the  electrons  NT p^}  times  the 

velocity  of  the  electrons  px/m  summed  over  those  values  of  momentum  from 

Px  to  infinity: 
o 

Jo  =  ;p  (e/m)  Px  N<Px>  dpx  (6) 

xo 

The  subscript  "o"  denotes  the  absence  of  external  fields,  and  J  is 

o 

referred  to  as  the  zero-field  current  density.  After  integration. 
Equation  6  becomes 

Jo  =  (4n  mek2T2/h3)  exp  T-e$/kT}  [A/m2]  (7) 

which  is  the  Richardson-Dushman  equation  (Reference  6). 
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Equation  7  is  more  commonly  expressed  as 

JQ  =  AT2  exp  ( -e<(>/kT )  (8) 

where 

A  =  4nmek2/h3  (9) 

and  A  has  a  numerical  value  of  approximately  120  [A/(cm2  °K2)]. 

If  the  electrons  experience  an  external  retarding  potential  of  height 

V$  above  f>  as  shown  in  Figure  1,  fewer  electrons  will  escape  so  the 

current  J  escaping  from  the  cathode  will  be  less  than  J  .  The  derivation 

o 

of  Equation  7  applies  with  the  use  of  ;■  +  Vs  instead  of  Then 
Equation  7  becomes 

J  =  (47rmek2  T2/h3)  exp  {-e  (<j>+Vs)/kT}  (10) 

where  the  zero  field  subscript  is  omitted.  Equation  10  can  also  be 
wri tten 

J  =  JQ  exp  i  -eVs/kT]  (11) 


where  JQ  is  given  by  Equation  7. 

Since  the  Fermi  level  Ep  is  measured  relative  to  the  bottom  of  the 
conduction  band,  current  will  flow  in  an  external  circuit  connected 
between  the  cathode  and  an  anode  unless  the  Fermi  levels  are  at  the  same 
potential.  So,  the  applied  voltage  (or  measured  voltage)  of  the 

a  pp 

anode  with  respect  to  the  cathode  is  just  the  difference  in  potential 
between  the  Fermi  levels  of  the  anode  and  cathode.  Figure  2  shows  for 
the  case  where  the  potential  minimum  occurs  at  a  collecting  body  (or  anode) 
that  the  applied  voltage  is 


(*c  ~  *a} 


The  derivation  through  Equation  11  applies  to  particles  of  charge  +e, 
although  electrons  are  conventionally  assigned  a  negative  charge.  This 
conflict  causes  no  trouble  if  the  direction  of  increasing  potential  is 
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taken  downward  in  Figures  1  and  2  for  the  development  through  Equation  11. 
However  when  measurements  are  involved  such  as  Vapp  which  assume  the 
negative  charge  for  the  electron,  the  direction  of  increasing  potential 
must  be  taken  upward.  For  this  reason  Vs  appears  with  a  negative  sign 
in  Equation  12. 

For  the  condition  that  the  potential  minimum  occurs  at  the  anode, 
combination  of  Equations  11  and  12  gives 

J  -  JQ  exp  {e(Vapp  +  *c  -  *a)/kT}  (13) 


which  shows  that  the  current  is  an  exponential  function  of  the  applied 

voltage.  This  is  known  as  the  "retarding-f ield"  or  "exponential"  range 

of  the  J-V  curve.  A  plot  of  1 n { J }  as  a  function  of  V  gives  a 

app  ^ 

straight  line  with  slope  e/kT.  The  current  J  equals  JQ  when 


(14) 


Figure  3  shows  a  schematic  of  the  I-V  Curve  under  these  conditions. 

Note  that  Equation  13  is  actually  independent  of  ,  as  long  as  J  is 
less  than  J  ,  since  d>  is  cancelled  by  the  J  term.  This  is  shown 
through  the  combination  of  Equations  8  and  13  which  yields 

J  =  AT2  exp  fe(Vapp  -  ^)/kT}  (15) 
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SECTION  II 
SPACE-CHARGE  THEORY 

Figure  1  shows  a  schematic  of  the  space-charge  region  outside  an 
emitting  surface.  The  position  and  magnitude  of  the  potential  minimum 
determines  the  amount  of  current  which  is  collected  at  an  anode  spaced 
some  distance  away  from  the  cathode.  If  the  potential  minimum  occurs  at 
the  anode,  the  potential  decreases  at  each  point  from  the  cathode  to  the 
anode.  If  the  potential  minimum  occurs  at  the  cathode,  the  potential 
increases  at  each  point  from  the  cathode  to  the  anode.  For  the  inter¬ 
mediate  case  of  a  potential  minimum  somewhere  between  the  cathode  and 
anode,  retarding-field  conditions  occur  near  the  cathode  and  accelerating- 
field  conditions  apply  near  the  anode.  For  this  intermediate  case  the 
electric  field  -dV/dx  at  the  minimum  is,  of  course,  zero.  Since,  for  the 
intermediate  condition,  the  current  is  greatly  influenced  by  the  mutual 
repulsion  of  the  electrons,  this  is  known  as  the  space-charge  region  of 
the  J { V I  curve. 

1.  DERIVATION  OF  LANGMUIR'S  DIFFERENTIAi  EQUATION 

The  following  discussion  derives  the  basic  equations  which  describe 
the  effect  of  space-charge  on  JiVi.  Electrons  with  energy  less  than 
Ep  +  e;  are  turned  back  right  at  the  metal-vacuum  interface.  Electrons 
with  energy  greater  than  Ep  +  e:  enter  the  space  outside  the  cathode  and 
are  acted  on  by  the  external  potential.  At  a  position  x  outside  the 
cathode  surface,  electrons  of  a  given  energy  from  the  momentum  distribution 
of  Equation  4  have  lost  momentum  due  to  the  barrier  at  the  surface  and 
due  to  the  potential  field  outside  the  surface.  As  indicated  in 
Figure  1,  which  shows  a  schematic  of  the  space  charge  region  outside  an 
emitting  surface,  the  net  loss  of  energy  is  Ep  +  e$  +  eVs  -  eVIxl  at 
position  x.  Note  that  Vlxt  is  measured  relative  to  the  potential 
minimum.  If  Equation  4  is  adjusted  for  this  lost  energy,  it  gives  the 
number  density  of  electrons  as  a  function  of  momentum  at  position  x  as 
written  in  Equation  16 

Nip  } dp  =  (4nmkT/h3)  exp  i Ep/kTlexp{-(p2/2m  +  Er  +  ef+eV  -eVIxl )/kT}dp 

X  X  r  X  r  j  X 

(16) 
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To  the  anode  side  of  the  potential  minimum,  the  lowest  energy 
electrons  at  position  x  will  have  energy  +  eVlx)  due  to  electric  field 
acceleration.  To  the  cathode  side  of  the  minimum,  electrons  will  be 
present  with  momentum  from  zero  to  infinity  toward  the  anode.  In 
addition,  electrons  will  be  present  with  energy  from  0  to  eVlx)  returning 
toward  the  cathode  from  the  direction  of  the  potential  minimum.  Then  to 
find  the  total  number  density  of  electrons  at  x  the  limits  of  integration 
will  be  +  (2meV<x')  '  with  the  +  sign  taken  in  the  anode  space  (to  the 
anode  side  of  the  potential  minimum)  and  the  -  sign  taken  in  the  cathode 
space.  This  integration  of  Equation  16  then  gives 

NT  V( x  1 1  =  [(?-nkT}3/2/h3]  exp  { -e( Vs+<f>)/kT fexpi eV{ x)/kT) 


•[1  +  Pi (eVlxl/kT)172}] 

(17) 

u  . t 2 

where 

Pi u }  =  (2//tT)  f  e  dt  =  error  function 

0  _ ** 

(18) 

integral 


In  Equation  17  the  -  sign  applies  in  the  anode  space  and  the  +  sign 
applies  in  the  cathode  space. 

The  electric  potential  outside  the  metal  is  governed  by  Poisson's 
equation 

d2V( x}/dx2  =  eN/e0  (19) 

where 

Q 

e0  =  1/(36tt  x  10  )[C/Vm]  =  permittivity  of  vacuum 

and  x  is  measured  from  the  position  of  the  potential  minimum.  The 
combination  of  Equations  17  and  19  and  the  introduction  of  the  normalized 
variables 

n  =  eVfxl/kT  (20) 

and 

f  =  x/x-j  (21) 
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where 

X,2  =  eQ  kTh3  exp  {e( Vs+ri> ) / kT > /2e2 ( 2-rrmkT ) 3/ 2  (22) 

yields  the  equation 

d2n/df.2  =  1/2  exp  {ni[l  7  P  (AT  }]  (23) 

This  equation  represents  the  potential  distribution  as  a  function  of 
distance  from  the  potential  minimum.  The  normalizing  constant  x^  as 
more  conveniently  written 

=  c0  (kT)3/2/Je  (M1/?  (24) 


where  J  is  given  by  Equation  10. 

Multiplication  of  Equation  23  by  dn/df;  and  integration  with 
respect  to  c  gives 

(dn/df. ) ^  -  n0^  =  [1  +  PiATl]  exp  In)  -  [1  +  2(nA)^2]  (25) 


where  h  is  defined  as 
o 


*0  H  dn/dc!r=o 


(26) 


As  before,  the  -  sign  of  the  +  applies  in  the  anode  space  and  the  + 
sign  applies  in  the  cathode  space.  In  order  to  maintain  a  positive 
direction  for  £,  from  cathode  to  anode,  the  square  root  of  Equation  25 
must  be  taken  with  a  negative  sign  in  the  cathode  space  and  positive 
sign  in  the  anode  space,  giving 

anode  dn./d 5  =  (h  2  +  [1  -  PiATl]  exp  fnJ-02(n/TT)1/2])1/2  (27) 

space  a  0  a  a  a 


1/2-H-1/2 


cathode  d£/dn„  -  -(n  +  [1  +  Pl^nT)]  exp  In  )~[l+2(n  / tt )  lfC~\) 

space  '  c  c 


(28) 
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The  subscripts  on  n  and  n  are  added  to  distinguish  the  anode  space  and 

3  C 

cathode  space  respectively.  The  reciprocal  expression  d;/dnc  is  used  in 
the  cathode  space  since  the  following  development  of  the  J-V  curve 
requires  n  as  a  function  of  ;  and  requires  ;  as  a  function  of  r 

a  C 

Rittner  (Reference  14)  has  compiled  a  collection  of  accurate  approximation 
formulas  for  the  solutions  to  Equations  27  and  28  for  A  =  0.  Rittner 
claims  a  maximum  of  .01%  relative  error  within  the  cited  range  for  the 
following  formulas: 

0  <  nc  1  0.20;  ;£  =  -  2nc1/2  +  .376126r,c  +  .0251956nc3/2 

+  . 000691 63n  2  -  .0009409n  5/2  (29) 

c  c 

-6. 1561 4 3r>  -65.9123n  2-59.4892n  3+4.103353r,  4 

0.20  <n  <3.0;  ;  =  - - - ^ 

-0.2629089  -14. 62018nc-49.59515r,cz-16. 06387  ,^+n^ 

(30) 


3.0  <  nc,  r,c  =  -2.55389  +  /~2  exp  {-r»c/2>  -  0.0123  exp  { -r>c} 

+  (l/3/T)[(nc/Tr)1/2  +  1]  exp  {-3nc/2)  (31) 

0  i  n  <  0.05;  n  =  0.25  ; }-  0.047016  ;  3  +  0.012627  ;  4  -0.003820;8  (32) 

a  a  a  a  d  3 

0.05  <  n  <  11.0;  n,  =  0.0003053192  -0.0010857509  ;  +  .2509922;  2 

a  a  a  a 

-0.04648354;  3  +0.01094699  ;  4  -0.002286296  £  5 

a  a  a 

+0.0003642352  E6  -0.00004062638  ;  7  (33) 

a  a 

+0.000002943653  ;  8-0. 000000123597  ;  9 

a  d 

+0.0000000022744  ;,10 

a 
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11.0  <  n  :  n  =  [('  +0.5088)/! ,2552]4/3  -1.5232(f  +0.5088)2/3 

—  a  a  a  a 

+  1.347  -0.495  (e  +0.5088)"2/3  -0. 140(£  +0.5088)"4/3 


Equation  31  comes  from  Langmuir  (Reference  7),  and  Equation  34  is  the 
inverse  of  a  series  given  by  Langmuir.  Equations  29  and  32  are  series 
solutions  to  the  differential  equations.  Equations  28  and  27  respectively, 
taken  from  Kleynen  (Reference  10).  Rittner  developed  Equations  30  and  33 
by  applying  standard  curve  fitting  techniques  to  tables  of  n  and  r  values 
given  by  Langmuir.  Langmuir  generated  the  n,  f,  values  by  numerical 
integration  of  the  differential  equations.  (Rittner  gives  a  different 

5 

coefficient  for  r  than  the  value  shown  here  from  Kleynen.  The  dis¬ 
cussions  in  '  ubsection  2  show  how  to  verify  that  Kleynen's  value  is  the 
correct  one. ) 

For  space-charge  conditions,  as  discussed  above,  the  potential 
minimum  occurs  between  the  cathode  and  anode  and  dV/dx  is  zero  at  the 
minimum.  For  this  reason  one  has 


dn/df’WV° 


for  space  charge  conditions.  A  plot  of  n(U  is  shown  in  Figure  4  along 
with  the  first  and  second  derivatives.  These  derivatives  will  be  useful 
below  in  discussions  of  the  behavior  of  JiVf.  Figure  4  was  generated 
through  the  use  of  Equations  29  through  34  and  Equations  23,  27,  and  28. 
Notice  that  f,  approaches  the  limiting  value  of  -2.55389  as  shown  by 
Equation  31  when  nc  tends  to  «. 

2.  SERIES  SOLUTIONS  FOR  THE  n,  ■  FUNCTIONS 

Although  Kleynen  (Reference  10)  does  not  discuss  how  he  derived  the 
series  solutions  for  the  n»  ■  functions  given  in  Equations  29  and  32,  the 
following  discussion  gives  one  way  of  deriving  the  terms  which  Kleynen 
gives  and  of  deriving  additional  terms: 

Let  h0  =  0  and  let 
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then  Equation  27  becomes 

du/d£  =  ([1  -  Ptui]  exp  {u2}  -  [l-2u//7"])1/2/2u  (37) 

From  Equation  18  the  derivative  of  P { u }  is 

dP/du  =  (2//T  )  exp  { -u2 }  (38) 

so  the  derivative  of  Equation  37  becomes 

d2u/dr,2  =  (  (u2-l)[l  -  Plu}]  exp  1  u2 }  +  [1  -2u/^T ] )/4u3  (39) 

Successive  derivatives  give 

d3u/d£3  =  ([2u4-3u2+3][l-P{u}]exp{u2)-[3-6u//V  +  2u3//_jT  ]) 

•(du/dc)/4u4  (40) 

and 

d4u/dr4  =  (d3u/dr,3)  (d2u/d;;2)/ (du/df, ) 

+([4u6-6u4+12u2-12][l-P{u)]exp{u2}-[4u5+8u3-24u]//7"  +12) 

•(du/ds)2/4u5  (41) 

Application  of  L'Hospital's  rule  (Reference  18)  to  (du/dc)2,  d2u/d£2 

3  3  4  4 

and  the  zero/zero  terms  in  d  u/df;  and  d  u/d£  gives  the  following  values 
for  these  derivatives  at  £  =  0: 

du/dc  =  1/2;  d2u/d?2  =  -1/6AT  ;  d3u/d£3  =  1/16 

d4u/dt4  =  -29/240  /T  (42) 

The  Taylor  series  (Reference  18)  about  £  =  0  then  gives 

utf.l  =  a1  e  +  a2  £2  +  a3  £3  +  a4  £4  (43) 
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where 

a]  =  1/2;  a2  =  -1/12/7  ;  a3  =  1/96;  a4  =  -29/5760.7  (44) 

The  coefficient  aQ  =  u{ t.  =  0}  is  zero  since  the  reference  voltage  and 
hence  u  is  chosen  to  be  zero  at  6  =  0. 

2 

To  obtain  an  equation  of  the  form  of  Equation  32  requires  n  =  u  , 

a 

that  is: 

na  =  b2  £2  +  b^  £3  +  b^  £4  +  b^  £3  (45) 

The  coefficients  bn  are  found  from  a  product  of  Equation  43  with  itself. 

2  2 

b^  =  a-|  ;  b3  =  2a^a^-,  b4  =  2a^a3  +  a^  ;  b&  =  23^  +  2a^a^  (46) 

Here  bQ  and  b-j  are  zero  since  aQ  is  zero.  Combination  of  Equations 
44  and  46  gives: 

b2  -  1/4;  b3  =  -1/12 77  ;  b4  =  (1/2  +  1/37/48  (47) 

b5  =  -39/5760v7“ 

Equation  32  results  from  a  substitution  of  the  numerical  values  from 
Equation  47  into  Equation  45. 

In  order  to  establish  Equation  29,  let 

v  =  77  and  £  =  -z  (48) 

Then  the  derivatives  of  v  with  respect  to  z  are  very  similar  to  those 
shown  in  Equations  37  through  41.  In  fact,  the  numerical  values  at 
z  =  0  of  dnv/d  zn  are  the  same  as  the  values  of  dnu/d£n  with  the 
variations  appearing  only  in  the  sign.  The  values  of  the  derivatives 
at  z  =  0  are  shown  in  Equation  49: 

dv/dz  =  1/2;  d2v/dz2  =  1/6/7  ;  d3v/dz3  =  1/16; 
d4v/dz4  =  29/240/7  ;  d5v/dz5  =  l/10n  +  9/128  (49) 
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The  Taylor  series  about  z  =  0  then  gives 

v(z!  =  CjZ  +  c2z2  +  c3z3  +  c4z4  +  c^z5  (50) 

whp  rp 

cn  =  ( 1 /n ; )  (dnv/dzn)  (51) 

The  function  vizi  is  then  obtained  by  substituting  £  =  -z  back  in 
Equation  50: 

vtr.l  =  g^  +  g2c2  +  g3£3  +  g4£4  +  g&£5  (52) 

where 

91  =  -1/2;  g2  =  1/12/T;  g3  =  -1/96;  g4  =  29/5760/7"  (53) 

g6  =  -(1/IOtt  +  9/1 28 )  / 1 20 

A  comparison  of  Equations  52  and  53  with  Equations  43  and  44  show  that 

ufrl  =  -  v  { r. }  (54) 

where  the  equality  is  taken  in  the  sense  of  the  Taylor  series  representation. 
Consequently  the  square  of  Equation  54  shows  that  the  Taylor  series  repre¬ 
sentation  for  n a { f, >  is  exactly  the  same  as  the  Taylor  series  representation 
of  n.lK)  at  r.  =  0.  This  shows  that  the  function  n,  taken  as  an  extension 
of  the  function  r>  is  continuous  at  £  =  0  and  has  continuous  derivatives. 

To  obtain  an  equation  of  the  form  of  Equation  29  requires  z T v >  rather 
than  vizi ,  that  is 

z  =  djv  +  d2v2  +  d3v3  +  d4v4  +  d5v5  (55) 

The  derivatives  of  the  inverse  function  z { v }  can  be  found  from  formulas 
relating  the  derivatives  of  the  inverse  function  to  the  derivatives  of 
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the  original  function,  some  of  which  are  given  by  Pierce  (Reference  19, 
p  100).  The  formulas  necessary  in  this  case  are 


dz/dv  =  1/v  (56) 

d2z/dv2  =  -v/(v)3  (57) 

d3z/dv3  =  [ 3 ( V ) 2  -  'v' v]/(v)5  (58) 

d4z/dv4  =  [10  v  v  V  -  15  (v)3  -  (v)2  v^v]/(v)3  (59) 

d3z/dv3  =  10[(  v  )2  (v)2  +  15  v  vlv  (v)2  -  105  (v)2  v  v 

+  105  (V)4  -  vV  (0)3]/(v)9  (60) 

where,  for  example,  d3v/dz3  has  been  shortened  to 'v,  and  d4v/dz4 
has  been  shortened  to  v1V.  Using  the  formula 

dn  =  ( 1/n! ) (dnz/ dvn )  (61) 

with  Equations  49  and  56  through  60  gives: 

d|  =  2  ;  d2  =  -  2/3.T"  ;  d3  =  -  (1  -  8/3tt )/6 ; 


d4  -  -2  ( 5/ 3tt  -  21/40)/9.>  ;  d5  --  (1/2  -72/5^  +  1  120/27t,2)/120  (62) 

Substitution  of  Equation  48  into  Equation  55  gives 


r  _  A  1/2 
*  ~  -  dlnc 


■  .3/2.2.  5/2 

'  d2\  '  3  ’c  ’  d4nc  '  d5nc 


(63) 


Finally,  the  use  of  Equation  62  in  Equation  63  gives  the  series  shown  in 
Equation  29. 


For  the  case  where  n  =  0,  Taylor  series  do  not  exist  for  n  If.)  or 

O  d 

£,  ■  d  1  •  This  is  seen  on  examination  of  the  third  derivative  of  Equation  27 
d  n  /  d  f 3  =  (1/2)[(1  -  P  { » 1 1  i)  exp  ini  -  T  /  v^rTn  ] 

o  add 

•[n2  +  (1  -  p  i,e  t)  expin  1  -(1-2 ,^)]1/2  (64) 

U  d  d  d 

The  term  1/  >rrT  increases  without  bound  as  n  approaches  zero,  showing 

d  d  J 

that  the  third  derivative  of  n  if.l  is  unbounded  at  f,  =  0.  Equations  23 
and  27  show  that  at  i  =  0 

dn  /df,  =  6  and  d2n,/d62  =  1/2  (65) 

d  U  d 
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The  calculation  of  the  derivatives  as  in  Equation  69  and  their 
evaluation  at  u  =  >  =0  gives 


v1  =  0; 


V2  =  17V 


»3  -  0;  v„  ■  -1/4^; 


(71) 


»5  '  4/1 5-  m  30, 


(2/n4  -3/.'.*)/24;  »?  "  2(4/^  -  1 5/^)/10b/7 


The  combination  of  Equations  67  and  6B  and  omission  of  the  zero  terms  gives 


v2na  +  v4  na  +  v5 


5/2  3  7/2 

1  +  V/-  n  +  V-7  n 

3  o  d  id 


(72) 


Figure  5  shows  a  plot  of  n.  versus  c  generated  from  Equation  72  with 

a 

nQ  =  0.5.  The  series  solution  with  up  to  the  seventh  power  of  u  from 
Equation  72  is  shown  along  with  the  series  solution  truncated  at  the 
sixth  power  of  u.  The  full  series  solution  is  useful  up  to  c  equal  to 
about  0.6.  Since  the  coefficients  involve  reciprocal  powers  of  n  in 
Equation  71,  one  might  expect  the  series  to  behave  better  at  larger 
values  of  nQ.  Although  the  series  is  able  to  handle  larger  values  of 
n.  as  6  becomes  large,  the  range  of  accuracy  for  c  is  still  1.  Note 
that  n  is  a  steeply  rising  function  of  r,  for  large  values  of  n  .  The 

3  O 

function  n  fr]  calculated  by  numerical  integration  is  shown  in  Figure  5 
for  comparison  with  the  series  solution.  The  numerical -integration  solution 
procedure  will  be  discussed  in  the  next  section. 

3.  NUMERICAL  SOLUTIONS  FOR  THE  n,  £  FUNCTIONS 

If  f,  =  0,  integration  of  Equation  28  for  tin  )  cannot  be  started 
at  n  =  f.  =  0  since  df/dnc  is  unbounded  at  this  point.  Most  numerical 
solution  techniques  for  differential  equations,  such  as  the  Runge-Kutta 
methods,  can  be  used,  however,  if  Equation  29  is  used  to  start  the  solution 
procedure.  The  function  njL  I  also  has  starting  problems  with  Runge- 

a 

Kutta  techniques. 

The  problem  with  n,(U  is  that  dn  /d£,  is  a  specified  function  of  n 

3  3  3 

starting  at  n  =0.  Any  solution  technique  such  as  the  Runge-Kutta  types 
which  extrapolate  values  of  n3  based  on  dn  /dc  will  always  give  n,td  =  0 

3  3  3 

for  all  c  if  it  begins  with  na(0}  =  0.  This  problem  is  solved  through  the 

3 
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use  of  a  Taylor  method  of  second  order  or  higher.  A  method  which  has 
been  very  useful  is  the  mixed  Runge-Kutta-Taylor  method  (Reference  20, 

3 

p.  275).  This  method  yields  convergence  of  order  h  ,  where  h  is  the  step 
size,  and  uses  the  first  and  second  derivatives  of  n  from  Equations  23 

d 

and  27.  This  solution  method  is  especially  convenient  since  the  second 

derivative  from  Equation  23  has  a  simple  form.  The  mixed  Runge-Kutta- 

Taylor  method  for  n  generates  a  sequence  of  numbers  n  ,  for  a  fixed 
3  n 

value  of  h,  by  the  recurrence  relation 

oQ  =  0  (73) 

nn+l  =  ”n  +  hG3tnn;h} 

wi  th 

G,u  ;  hi  =  fir,  }  +  (h/2)  f{n  +  hf{n  }/3)  (74) 

o  3  3  da 

where 

f  =  dn  /dr  (75) 

a 

and  f  and  f  =  df/dr,  are  given  by  Equations  23  and  27. 

The  accuracy  of  the  Runge-Kutta-Taylor  method  solving  for  n Jr,}  is 

a 

limited  less  by  the  step  size  h  than  it  is  by  the  accuracy  of  the  values 
calculated  for  the  error  function  Piul.  The  most  satisfactory  representation 
for  P { u 1  used  in  this  work  comes  from  Abramowitz  (Reference  21,  page  299, 
Equation  7.1.26).  This  representation  is 

P{ ul  =  1  -  (a^t  +  +  a^t^  +  a^t^  +  a^t^)  exp  {-u^l  (76) 

where 

t  =  1/(1  +  pu)  (77) 

and 

p  =  .3275911  a1  =  .254829592 

a2  =  -.284496736  a3  =  1.421413741  (78) 

a4  =  -1.453152027  a5  =  1.061405429 
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Since  the  machine  available  for  calculations  in  this  work  had  only  a  six 
digit  capability,  the  numerical  values  in  Equation  78  were  rounded  to 


p  =  .327591 

a3  =  1.42141 


a1  .25483  a2  =  -.2845 

a4  =  -1.45315  a5  =  1 .06141 


(79) 


It  is  important  in  this  round-off  procedure  to  have  the  terms  through 
a^  to  add  to  exactly  1.  For  example  the  rounding  of  a2  to  -.284497  instead 
of  the  value  in  Equation  79  results  in  greatly  increased  error  in  the  n 

d 

values  near  e,  =  0.  The  use  of  Equation  76  for  P { u >  along  with  the 
numerical  values  in  Equation  79  in  the  Runge-Kutta-Taylor  solution 
procedure  generates  accurately  the  values  of  n  { e, >  tabulated  by  Kleynen. 
With  this  representation  of  P{ul,  step  size  h  less  than  h  =  .01  does  not 
increase  the  accuracy  of  the  solution.  Step  size  h  =  .01  generates  the 
Kleynen  (Reference  10)  tables  for  n  { r, }  accurate  to  the  fourth  place 
for  f,  ^  5,  and  h  =  1.0  gives  good  accuracy  above  K  =  100.  A  BASIC 
computer  language  program  using  Equations  73  through  79  for  solving 
Equation  27  is  shown  in  Figure  6.  A  sample  print-out  from  the  program 
for  h  =  .01  is  shown  in  Figure  7. 

For  the  case  where  f,  is  greater  than  zero  a  first  examination  seems 
to  indicate  that  a  Runge-Kutta  method  would  be  better  for  solving  for 
n a { }  than  the  Runge-Kutta-Taylor  method  of  Equation  74.  Since  the  third 
derivative.  Equation  64,  of  o  l£}  is  unbounded  at  f,  =  0,  one  might 

a 

expect  problems  with  the  calculation  of  »  from  Equations  73  and  74  which 

a 

use  values  of  the  second  derivative  of  n  .  Note  that 

d 

fin  1  =  dZn  Jdr2  (80) 

a  a 


This  concern  is  based  on  the  Taylor  Theorem.  When  a  function  is 
generated  from  its  first  and  second  derivatives,  the  remainder  term  is 
related  to  the  third  derivative.  In  addition,  when  n  is  greater  than 
zero,  the  Runge-Kutta  methods  do  not  generate  the  spurious  solution 
n.  =  0  as  discussed  earlier. 

d 
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However,  in  practice,  the  Runge-Kutta-Taylor  method  yields  much 
quicker  convergence  than  the  simplified  Runge-Kutta  method  which  uses 
the  same  number  of  function  evaluations  per  step.  This  superior  conver¬ 
gence  is  obtained  over  a  wide  range  of  values  of  n  .  The  simplified 
Runge-Kutta  method  is  described  by 


n 


o 


=  0 


wi  th 


Vl  =  nn  +  hK2{Vh} 


K_in  ;h>  =  (l/2)[f{ n_  >  +  flna  +  hfln  M] 

c  a  a  a  a 


(81) 

(82) 


where 

f  =  dn  /dr,  (83) 

a 

and  f  is  given  by  Equation  27.  The  convergence  of  the  Runge-Kutta- 

Taylor  method  compared  to  the  convergence  of  the  simplified  Runge-Kutta 

method  is  shown  in  Table  1  for  the  example  hQ  =  .003  at  5  =  5.0.  Later 

work  will  show  that  these  conditions  correspond  approximately  to  a 

cathode  with  I  30  mA/crn^,  T  =  650°C,  I  =  1 . 002 1 o ,  and  cathode-anode 

spacing  0  =  .05  mm.  Table  1  also  compares  the  convergence  of  these  two 

methods  for  the  example  n  =  .5  at  £  =  1000.  Later  work  will  show  that 

o  c  2 

these  conditions  correspond  approximately  to  a  cathode  with  I  =  5A/cm  , 

T  =  900°C,  1=1.11,  and  0  =  1.0  mm. 

o 

The  superior  convergence  of  the  mixed  Runge-Kutta-Taylor  method 

3 

results  from  its  third  order  convergence,  order  of  h  ,  compared  to  the 

2 

second  order  convergence,  order  of  h  ,  for  the  simplified  Runge-Kutta 
method.  For  small  values  of  the  simplified  Runge-Kutta  method 
generates  tne  solution  n  =  n  r  which  has  nearly  the  error  of  the  spurious 
solution  n  =  0  generated  for  the  case  h  =0.  On  the  other  hand,  the 

3  O 

third  derivative  of  n  is  large  over  a  very  small  interval  of  values  of  £ 
for  small  h  .  For  this  reason,  the  total  change  in  the  second  derivative 
is  not  large  near  £  =  0  and  the  Runge-Kutta-Taylor  method  accurately 
calculates  na •  For  large  values  of  h0  the  value  of  na  is  dominated  by  the 
first  derivative  so  that  large  variations  in  the  second  derivative 
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TABLE  1 

COMPARISON  OF  THE  CONVERGENCE  OF  THE  MIXED  RUNGE-KUTTA-TAYLOR ,  R-K-T, 


METHOD  WITH  THE 

SIMPLIFIED  RUNGE- 

KUTTA,  R-K,  METHOD 

nQ  =  1  x  10-3  £  =  5 

step  size 

n  from  R-K-T 
a 

n  from  R-K 
a 

.005 

3.61135 

3.61046 

.01 

3.61139 

3.60858 

.05 

3.61163 

3.58151 

.1 

3.61215 

3.53811 

.5 

3.63247 

3.13272 

1.0 

3.69794 

2.62286 

5.0 

5.79598 

.300455 

nQ  =  0.5  £  =  1000. 

step  size 

n  from  R-K-T 
a 

from  R-K 
a 

.5 

7287.76 

7286.99 

1 

7288.19 

7286.38 

5 

7291 .94 

7276.32 

10 

7297.70 

7259.18 

100 

7494.38 

6838.45 

200 

7804.92 

6328.12 
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near  5  =  0  do  not  cause  problems  with  the  accuracy.  Consequently,  the 
mixed  Runge-Kutta-Taylor  method  accurately  calculates  na  also  for  large 
h  .  Since  the  mixed  Runge-Kutta-Taylor  method  converges  faster  for 
small  nQ  and  at  least  as  fast  for  large  hQ  it  is  used  for  the  calculations 
reported  later  in  this  work.  The  criterion  for  the  choice  of  h,  used  for 
nQ>0  is  the  same  as  the  criterion  for  hQ  -  0  discussed  above. 

4.  THE  CURRENT  -  VOLTAGE  RELATIONSHIP  FOR  J  <_  JQ 

If  J  <  J  ,  Equations  27  and  28  give  the  voltage  according  to  the 
procedure  outlined  by  Langmuir  (References  7  and  8).  If  the  space- 
charge  potential  minimum  occurs  at  the  anode,  the  voltage  can  be 
calculated  from  Equation  13.  Otherwise,  the  potential  minimum  occurs 
somewhere  between  the  cathode  and  the  anode,  and  the  value  of  the 
potential  minimum  relative  to  the  metal  surface  is  V$  from  Equation  11: 

V$  =  (kT/e)ln  (84) 

The  distance  from  the  potential  minimum  to  the  cathode  is  found  from 
Equation  28  since  at  the  cathode  nc  =  (e/kT)Vs: 

£  =  £  in  )  =  £rUn{J  /J}}  (85) 

La  L  U 

cathode 

Recall  that  £c  and  nc  are  measured  relative  to  the  potential  minimum  and 
x  =  £  x^  from  Equation  21.  The  distance  from  the  potential  minimum  to 
the  anode  is  then  the  anode  -  cathode  spacing  D  minus  the  distance  from 
the  potential  minimum  to  the  cathode  given  by  Equation  85.  In  normalized 
terms,  then,  the  locations  of  the  anode  is  £,  where 

Q 

£  =  0/T/  x„  +  £ r{1n{J  /J}}  (86) 

a  CL  O 

Here  the  normalizing  constant  x^,  based  on  Equation  24,  given  by 

x22  =  eQ(kT)3/2  /e(2*m)1/Z  =  x^J  (87) 
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allows  the  current  density  J  to  appear  explicitly  in  Equation  86.  Then 
from  Equation  27  the  potential  Va  of  the  anode  relative  to  the  potential 
minimum  is 

Va  =  (kT/e)na{DAT  /x2  +  cctlntJQ/J}l  (88) 

Figure  8  shows  a  schematic  of  the  potential  distribution  between 
the  anode  and  cathode  under  space-charge-1 imi ted  conditions.  This 
diagram  shows  that  the  anode  to  cathode  voltage  V  is 

V  =  V  -  V.  (89) 

d  b 


and  the  applied  voltage  which  is  offset  by  the  Fermi  levels  is 


V 

app 


V  +  *a 


V 


C 


(90) 


(A  similar  situation  is  discussed  above  in  the  derivation  of  Equation  12.) 
Then  from  Equations  84,  88,  and  89 

V  =  (kT/e ) [ln{ J/JQ}  +  na(D /T  /x2  +  f,c(  ln{  JQ/J}}  ]  (91) 

where  hQ  =  0  in  Equations  27  and  28  and  D  is  the  cathode-anode  spacing. 

Equation  91  is  instructive  as  a  summary  of  the  solution  procedure 
outlined  by  Langmuir.  Each  term  in  Equation  91  represents  a  part  of  the 
physical  behavior  of  the  emission.  The  term  ln{J/JQ}  is  negative  and 
represents  a  voltage  offset  due  to  the  existence  of  the  potential 
minimum  V  .  The  term  r  { 1 n { J  /JH  is  negative  and  represents  the 
distance  from  the  potential  minimum  to  the  cathode.  The  term  Dv'T”  /x^, 
serving  as  the  argument  for  the  function  ri  ,  is  responsible  for  the 
basic  3/2  power  behavior  of  the  current  as  a  function  of  voltage. 
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5.  BEHAVIOR  OF  J{ V}  FOR  SPACE-CHARGE-LIMITED  CONDITIONS 

Three  characteristics  of  the  J{ V}  curve  are  of  special  interest  for 

space-charge-limited  conditions.  Each  of  these  three  characteristics, 

the  transition  from  retarding-f ield  to  space-charge-limited  operation, 

the  inflection  point,  and  the  transition  from  space-charge-limited 

operation  to  the  zero-field  condition  J  =  J  ,  can  be  used  in  an 

experimental  determination  of  the  zero  field  current  density  J  which  is 

o 

characteristic  of  the  cathode  under  test.  The  derivative  of  Equation  91 


dV/dJ  =  ( kT/eJ ) [1  +  (D/3~/2x2  -  d£/dnc 


)  dn/df; 
a 


(92) 


and  the  second  derivative 


d2V/dJ2  =  (kT/eJ2) 


-1  -  (dn  /dr, 

a 


)  [-dc/dn. 


a 


+  ST D/2x9] 

I  C. 


?  ? 
+  (d  n  /df, 

a 


f,1 

a 


)  [-df,/dnf 


+  ST  D/2 x2Y 


(93) 


+  (dn  / df, 

d 

)  [  -  (dc/dn 

r  1  L 

)3  (d2n  /dr2 

1  v. 

)  +  D/J“/4x2] 

r 

a 

nc 

nc 

show  more  explicitly  the  behavior  of  the  current-voltage  curve  for 
space-charge-limited  conditions.  Here  n'  and  £ 1  are  defined  as 

C  3 


and 


=  ln{ JQ/J} 

(94) 

=  ST  D/x2  +  Lc(ln{J0/J}} 

(95) 

2  2 

The  derivatives  dJ/dV  and  d  J/dV  are  obtained  from  Equations  92  and  93 

through  the  use  of  Equations  56  and  57.  Equation  57  is  also  used  to 

2  2 

express  the  derivative  d  f./dn  in  Equation  93. 
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Figure  9  shows  J{V}/JQ,  ( dJ / dV ) / JQ  and  (d^J/dV^)/JQ  for  an  example 
cathode  having  JQ  =  10  mA/cm^,  T  =  700°C,  and  D  -  0.2  mm.  The  curves  are 
normalized  by  the  factor  JQ  for  ease  of  plotting.  The  curves  in  Figure  9 
are  shown  dashed  in  the  retarding  field  range  and  solid  in  the  space-charge 
range.  The  retarding-f iel d  range  is  described  by  Equation  15. 

At  the  transition  from  the  retarding-field  range  to  the  space-charge 

range.  Figure  9  indicates  that  J { V }  is  a  smooth  curve  and  that  dJ/dV  is 

2  2 

continuous.  The  second  derivative  d  J/dV  ,  however,  is  discontinuous  at 
this  transition  point.  The  second  derivative  is  much  larger  on  the 
retarding-field  side  than  it  is  on  the  space-charge  side.  (In  the 
generation  of  Figure  9,  Equations  92  and  93  will  generate  very  curious, 
though  spurious,  values  when  applied  in  the  retarding-field  range.) 

Toward  increased  voltage  in  the  space-charge  range,  notice  that 
2  2 

d  J/dV  passes  through  zero,  indicating  an  inflection  point  in  the  J f V } 

curve.  Crowell  discusses  the  existence  of  this  inflection  point  through 

plots  of  the  results  of  many  numerical  calculations.  However  he  does  not 

relate  the  behavior  of  the  current- vol tage  curve  back  to  the  equation  for 

J{V}.  Crowell  also  does  not  point  out  that  for  practical,  high-current- 

density  cathodes  this  inflection  point  occurs  virtually  at  J  =  JQ.  Note 

that  the  example  in  Figure  9  is  for  a  cathode  with  very  low  current 
2 

density,  10  mA/cm  ,  and  a  very  close  diode  spacing  of  0.2  mm.  A  detailed 
discussion  of  this  inflection  point  follows  below  in  Subsection  b. 

Perhaps  more  important  than  the  existence  of  the  inflection  point  is 

the  fact  that  J { V }  tends  to  JQ  with  zero  slope.  Work  discussed  below  in 

Subsection  a  will  show  that  this  behavior  of  J{V)  near  Jq  is  characteristic 

of  high-current-density  cathodes  as  well  as  the  low-current-densi ty 

cathode  of  Figure  9.  The  work  below  and  the  discussions  in  Section  IV 

on  the  behavior  of  J { V }  in  the  accelerating-field  range  predict 

2  2 

discontinuities  in  dJ/dV  and  d  J/dV  at  J  -  JQ.  This  behavior  should 
show  up  as  rapid  variations  in  these  derivatives  for  cases  where  experi¬ 
mental  conditions  (such  as  edge  effects  or  emission  patches)  obscure 
details  of  the  J { V >  curve.  These  rapid  variations  in  the  derivative 
should  serve  as  a  better  indicator  of  the  value  of  JQ  than  the  inflection 
point  measurement  discussed  by  Crowell  (Reference  12). 
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a.  Transition  to  J  =  JQ 

At  the  point  where  J  is  just  equal  to  JQ  from  the  space-charge 
range,  the  space-charge  minimum  occurs  right  at  the  cathode.  In  this  case 
n ’  from  Equation  94,  which  represents  the  voltage  offset  due  to  the 
potential  minimum,  is  zero.  Then  in  Equation  92  since  dnc/dc  is  zero  at 
this  point,  the  term  dt/dnc  is  unbounded.  Then  dV/dJ  is  unbounded  which 
implies  that 

dJ/dV  =  0  at  J  =  JQ  (96) 

This  is  true  for  any  value  of  JQ.  This  behavior  appears  clearly  in  an 

example  treated  by  van  der  Ziel  (Reference  9).  (See  van  der  Ziel's 

2 

Figure  1  with  I$  =  0.01  A/cm  .) 

For  many  practical  problems  this  behavior  of  the  space-charge 
equation  only  applies  over  values  of  J  extremely  close  to  JQ  and  is  not 
readily  visible  in  the  current-voltage  plots.  For  example,  the  zero  slope 
of  J  as  a  function  of  V  at  J  =  JQ  is  not  visible  in  Figure  24  which 
treats  a  high-current-density  cathode.  However,  later  work  in  Section  IV 
shows  that  dJ/dV  is  greater  than  zero  at  J  =  JQ  from  the  accelerating- 
field  range.  This  discontinuity  in  the  first  derivative  and  the  accompany¬ 
ing  discontinuity  in  the  second  derivative  should  show  up  in  plots  of 
experimentally  taken  derivatives  even  when  the  transition  at  J  =  JQ  is 
masked  in  the  J { V }  plot  itself. 

b.  The  Inflection  Point 

The  inflection-point  current  is  defined  as  an  approximate 
measure  of  JQ  in  the  1950  I.R.E.  Standards  on  Electron  Tubes  (Reference  16). 
This  standard  is  based  on  the  observation  that  J{ V}  has  an  increasing 
slope  for  small  J  and  a  decreasing  slope  for  large  J.  The  inflection 
point  is  then  the  transition  point  from  increasing  to  decreasing  slope  of 
J { V > ,  the  point  where  d  J/dV  is  zero.  Crowell  (Reference  12)  was  the 
first  to  point  out  that  Langmuir  space-charge  theory  predicts  an  inflection 
point  for  some  J  less  than  JQ.  Crowell  used  a  numerical  study  for  many 
values  of  JQ  to  support  the  use  of  the  inflection  point  as  an  approximate 
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measure  of  JQ.  This  section  goes  beyond  the  work  of  Crowell  by  showing 
on  a  theoretical  basis  how  Langmuir  space-charge  theory  always  produces 
an  inflection  point  at  some  J  less  than  JQ. 

Equation  93  along  with  Equations  92  and  57  give  a  fairly 

2  2 

complicated  expression  for  d  J/dV  .  However,  in  the  limit  as  J  approaches 
J  ,  the  term  df./dn  is  unbounded.  Use  of  this  fact  allows  a  simple 
expression  for  d^J/dv  at  J  =  JQ 


2  2, 

d  J/dV 

|  J 


=  -(J/2)(e/kT)2(dr,a/d;  |  ) 


(97) 


At  J  =  Jfl,  f,  is  zero  so  from  Equation  95 

f '  -  /J”  D/x0 

a  c 


(98) 


Then  Equation  97  shows  that  the  second  derivative  of  J{V)  is  always 

2  2 

negative  and  finite  at  J  .  Although  Figure  9  indicates  that  d  J/dV  has 
a  large  negative  value  for  that  example.  Equation  97  shows  that  the 
second  derivative  is  finite. 

For  small  J,  in  the  retarding-f ield  range.  Equation  13  shows 
that  J { V }  and  its  derivatives  are 

J  =  JQ  exp  {e  (Vapp  +  <frc  -  <}>a)/kT}  (99) 

dJ/dV  =  (eJo/kT)  exp  le  (Vapp  +  ;c  -  Da)/kT}  (100) 

and 

d2J/dV2  =  (e/kT)2  Jo  exp  (e  (Vapp  +  ^  -  *a)/kT)  (101) 

Equation  101  shows  that  the  second  derivative  is  always  positive  for 

2  2 

small  J.  Comparison  of  Equations  97  and  101  shows  that  d  J/dV  changes 

sign  for  some  value  of  J  between  J  =  0  and  J  =  J  .  In  other  words, 

o 

there  is  an  inflection  point  of  J f V }  at  this  value  of  J. 

This  inflection  point  has  the  usual  properties  except  for  the 
case  of  very  small  values  of  v^J-  O/x^  discussed  in  Subsection  c.  In  this 
case  a  discontinuity  in  the  second  derivative  coincides  with  the  inflection 
point.  The  second  derivative  then  passes  discontinuously  from  positive 
values  to  negative  values. 
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c.  Transition  from  Retarding-field  to  Space-Charge-Limited 

Condi tions 

This  section  explores  the  detailed  behavior  of  Langmuir  space- 
charge  theory  at  the  transition  from  the  retarding-field  range  into  the 
space-charge-limited  range.  Van  der  Ziel  (Reference  9)  mentions  this 
transition  point  and  gives  it  the  name  "exponential  point".  However 
van  der  Ziel  does  not  discuss  the  behavior  of  the  derivatives  of  J { V } 
at  the  exponential  point,  probably  since  he  did  not  have  Equation  91  or 
an  equivalent.  The  examination  below  of  J { V]  and  its  derivatives  show 
that  the  J { V >  curve  is  smooth  at  the  transition  from  the  retarding  field 
region  into  the  space  charge  region  ("exponential  point").  In  addition, 
the  derivative  dJ/dV  is  continuous  but  not  smooth,  and  the  second  derivative 
is  discontinuous  at  the  exponential  point. 

For  the  condition  that  the  space-charge  just  begins  to  reduce  the 
current  from  the  value  it  would  have  in  the  retarding-field  range,  the 
space  charge  minimum  occurs  right  at  the  anode.  In  this  case  the  dis¬ 
tance  from  the  potential  minimum  to  the  anode  is  zero  so  the  value  of  £ 
at  the  anode,  designated  r'  above,  is  zero.  From  the  definition  of  £' 
in  Equation  95,  it  follows  that 


ccnn(Jo/Jen  -  -  ^  d/x2 


(102) 


where  J  is  the  current  density  at  the  exponential  point.  Then  since  n  , 

c  d 

evaluated  at  r  =  0,  is  zero,  Equation  91  reduces  to 


Ve  =  (kT/e)  In  {Je/JQ} 


(103) 


where  Vg  is  the  voltage  at  the  exponential  point.  Since  h  ,  which  is 


dn/d£  evaluated  at  £  =  0,  is  zero,  Equation  92  reduces  to 


dV/dJ! 


J  =  J 


=  kT/e  J, 


(104) 


and  Equation  93  reduces  to 

d2 V/dJ2 i  =  (kT/e  J  ' 2)[  -1  +  1/2  (-d£/dn.  +  /J~D/2x9)2]  (105) 

J  =  J  6  C  n' 

1  e  c 

2  2 

In  Equation  105  use  has  been  made  of  Equation  23  to  evaluate  d  n  /d£  at 

3 

^  =  °- 
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If  Equations  103  through  105  are  written  in  terms  of  J{ V j  instead 
of  V{ J }  they  become 

Je  =  JQ  exp  { e V  / k T }  (106) 

dJ/dV  -  ej  / k T  (107) 

•J  =  Je  " 
and 

d2 J/dV2  =  (e/kT)2  J  [1  -  1/2  (-d  ./&.  ,  +  .TO /2x?)2]  (108) 

'j  =  J  e  ^ 

e  1  c 

Comparison  of  Equations  99  and  100  with  Equations  106  and  107  shows 
that  J { V }  and  dJ/dV  are  both  continuous  at  the  exponential  point  and 
consequently  J { V }  is  smooth.  Note  that  V  and  V  -  :•  are  the 

d  pp  C  a 

same  from  Equation  90.  Comparison  of  Equations  101  and  108  shows  that 
2  2 

d  J/dV  is  always  greater  on  the  retarding  field  side  of  the  exponential 
point  than  it  is  on  the  space-charge  side.  From  Equation  28  the  term 
d.-'./dn  in  Equation  108  is  always  negative.  (This  is  also  seen  in 
Figure  4.)  Then  the  term  in  brackets  in  Equation  108  is  always  less 
than  the  value  calculated  from  Equation  101.  Then  the  second  derivative 
is  discontinuous  so  that  the  first  derivative  is  not  smooth  at  the 
exponential  point.  Figure  10  shows  a  detailed  sketch  of  dJ/dV  at  the 
exponential  point  for  the  example  of  Figure  9. 

For  many  cathodes  where  Jo  is  relatively  large,  ■'  cl  ln{  JQ/J  1 1 
is  very  nearly  its  minimum  value  of  -2.5539  before  it  is  equal  in 
magnitude  to  D»"J  /x^.  In  this  case  Jg  at  the  exponential  point  is 
much  less  than  JQ  and,  from  Equation  102,  Jg  is  given  by 

Je  *  (2.5539  x2/D)2  (109) 

In  this  case  of  large  JQ  then,  at  the  exponential  point  becomes 
nearly  independent  of  the  zero  field  current  JQ.  In  this  case  dt/dn 
at  is  very  small,  so  the  second  derivative  is  approximately 

d2J/dV2:  --  (e/kT)2  J  [1  -  (l/2)(?.5/2)2]  (110) 

1  1  -  1  “ 


32 


AFWAL -TR-81 -1 01 9 


which  is  approximately  a  factor  of  5  jump  in  the  second  derivative  at  the 
exponential  point.  For  the  example  of  Figure  9  the  second  derivative 
jumps  by  a  factor  of  5.7.  Note  in  Figure  9  that  the  second  derivative 
is  off  scale  to  the  retarding  field  side  of  the  transition  point. 

For  cathodes  with  very  small  JQ  the  term  AT  D/x^  in  Equation  102 
will  be  very  small  leading  to  a  very  small  K ' •  In  this  case  -dt/dr 

c  2  2  ^ 

would  be  quite  large  in  Equation  105  rendering  a  3/6 v  negative.  If 

2  2 

d  J/dV  is  negative  to  the  space-charge-1 imi ted  side  of  the  exponential 
point,  there  is  no  additional  inflection  point  in  the  J { V }  curve  and  the 

slope  of  J { V }  always  decreases  in  the  space-charge-1 imi ted  range.  Since 

2  2 

d  J/dV  is  positive  in  the  retarding  field  range,  this  discontinuity  to 

2  2 

negative  d  J/dV  can  be  quite  pronounced  for  low-current-density  cathodes. 


d.  Experimental  Use  of  the  Exponential  Point 

The  discontinuity  in  the  second  derivative  at  the  exponential 

point  should  lead  to  a  sharp  peak  in  the  second  derivative  of  experimental 

J{ V)  curves.  The  current  J  and  the  applied  voltage  1/  „  can  then  be 

e  r  app 

measured  at  this  peak  in  the  second  derivative.  Then,  if  JQ  is  not  too 

large,  the  theory  developed  in  Subsection  c  and  expanded  below  calculates 

the  zero  field  current  J  and  the  work  function  difference  4>  -  4>  - 

o  a  c 

A  suggested  limit  on  JQ  for  accurate  calculations  is  that  be  greater 
than  approximately  JQ/50.  This  condition  on  JQ  is  satisfied  for  large 
tj>  cathodes  or  by  sufficiently  low  temperature  T. 


Based  on  one  measurement  of  current  J  and  voltage  V  ,  the  zero- 

e  e 

field  current  JQ  of  the  cathode  and  the  difference  in  cathode  and  anode 

work  function  can  be  calculated  without  reference  to  the  Richardson 

Equation  8.  Equation  8  can  then  be  used  to  calculate  the  cathode  work 

function  from  the  zero-field  current  density.  From  Equation  90  the 

measured  voltage  V,  „  is  related  to  the  voltage  V„  by 
app  3  e  J 


V  =  V 
app  e 


(111) 


From  Equation  102  the  current  Jg  is  related  to  JQ  by 

Jo  =  Je  exp{nc{,/Je  d/x2}1 


(112) 
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and  Equation  103  relates  V  to  J  and  J  .  Then  Equations  103  and  112  can 

be  used  to  generate  plots  of  JQ  and  Vg  as  functions  of  Jg.  For  a  measured 

value  of  J&  and  with  measured  values  of  T  and  D,  the  zero  field  current 

density  can  be  calculated  from  Equation  112.  Then  Equation  103  can  be 

used  to  calculate  V  .  From  the  calculated  value  of  V  and  the  measured 
e  e 

value  of  V  ,  Equation  111  then  produces  the  difference  between  the 
app  r 

cathode  and  anode  work  functions  U  -  ). 

C  3 

The  problem  with  the  use  of  Equation  112  is  that  it  requires 
nc f t )  in  the  cathode  space.  A  method  for  calculating  this  function  could 
be  derived  in  a  straightforward  manner.  However,  since  the  inverse 
function  }  can  be  calculated  from  Equations  29  through  31,  it  seems 
easier  at  this  point  to  restate  the  problem,  so  that  this  inverse  function 
can  be  used.  The  use  of  this  inverse  function  is  possible  through  the 
construction  of  an  auxiliary  function  f { J }  where  based  on  Equation  102 

f{ J}  =  £c(ln{J0/J}}  +  0/T  /x2  (113) 

Then  finding  the  solution  J  to  Equation  112  for  a  given  JQ  is  equivalent 
to  finding  the  solution  to  Equation  113  such  that 

f{Je)  =  0  (114) 

This  allows  a  value  of  Jg  to  be  generated  for  each  JQ  over  a  range  of 
values  of  JQ.  This  can,  of  course,  still  be  plotted  as  JQ { }  for  use 
in  interpreting  experimental  data. 

A  plot  generated  from  the  BASIC  program  discussed  below  is  shown 
in  Figure  11  for  two  examples.  In  one  example  the  temperature  is  700°C 
and  the  spacing  is  0.2  mu  while  in  the  other  example  the  temperature  is 
900°C  and  the  spacing  is  1.0  mm.  Actually,  the  values  of  temperature  and 
spacing  are  not  unique  to  each  curve  in  Figure  11.  Rather,  the  family  of 
curves  suggested  in  Figure  11  depends  on  the  parameter  t  where 

o  z  A  1/2 

t  =  D/x2  =  (D/T3/4) (e/2irm  /cQ  k3/2)  (115) 

as  seen  in  Equation  113.  Figure  11  also  indicates  the  asymptotic  values 
where  Jg  ceases  to  increase  for  increasing  values  of  J  .  The  asymptotic 
value  is,  of  course,  dependent  on  cathode  temperature  and  diode  spacing 
and  is  given  by  Equation  109. 
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Equation  114  can  be  solved  for  Jg  by  an  iterative  technique 
known  as  Newton’s  method  (Reference  20,  p  277).  Based  on  an  initial 
guess  J-j  at  the  solution,  this  method  generates  a  sequence  of  approximations 
by  the  recurrence  relation 

Jn  +  1  =  Jn  -  f{Jn)/f{Jn}  n  =  1,2,3 .  (116) 

The  iteration  is  terminated  when  J  ,  is  sufficiently  close  to  J  and 

J  .  .  is  taken  as  the  value  of  J  .  Since  J  is  always  less  than 
n  +  1  ^  e  e  J 

(2.554  x^/D;  from  Equation  109,  this  is  a  convenient  starting  value  for 

J,.  However,  for  small  values  of  J  this  formula  gives  a  value  of  J 
1  o  e 

which  is  greater  than  JQ.  In  this  case  a  smaller  starting  value  for  J-j 
is  necessary,  for  example  J-j  =  0.5  J  .  In  Equation  116  the  notation 
f'  {J^l  represents  df/dJ  evaluated  at  J  ,  that  is 

f'{Jn)  =  (-l/J)(df./dnc)  +  D/2  x2  ,/3T  (117) 

where  df/dn  is  calculated  from  Equation  28  with  n  evaluated  at 

nc  =  ln{ JQ/Jn>  (118) 

A  BASIC  program  which  solves  Equation  114  and  generates  and  Vg  for  a 
range  of  values  of  JQ  is  shown  in  Figure  12. 

Another  way  of  representing  the  behavior  of  J { V }  at  the 

exponential  point  is  shown  in  Figure  13.  This  figure  shows  both 

and  Je  plotted  as  functions  of  V  .  The  calculations  were  made  by  the 

BASIC  program  of  Figure  12.  Figure  13  is  especially  convenient  for 

interpreting  measured  data  since  it  gives  both  JQ  and  V0  based  on  the 

measured  value  of  J  .  From  V  and  the  value  of  V  „  at  the  transition 
e  e  app 

point.  Equation  111  then  gives  the  difference  in  work  function  (4>  -  <j>  ) . 
Figure  13  shows  the  asymptotic  value  of  Jg  for  the  example  T  =  700°C 
and  D  =  0.2  nun.  Figure  13  also  shows  the  asymptotic  line  for  JQ(Ve} 
which  passes  through  the  asymptotic  value  of  Jg  at  Vp  =  0  with  a  slope 
of  kT/e.  The  existence  of  an  asymptotic  line  for  constant  is  seen 
from  Equation  103.  For  experimental  interpretation  families  of  curves  of 
JQ  or  Jg  as  functions  of  eV/kT  can  be  generated  just  as  the  sample  curve 
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was  generated  for  Figure  13.  With  the  V  axis  normalized  to  eV/kT,  the 
family  of  curves  is  generated  with  the  single  parameter  t  =  D/x^  just  as 
in  the  generation  of  the  family  of  curves  like  those  in  Figure  11. 

Although  the  representation  is  not  as  intuitive  as  the  families 

of  curves  suggested  in  Figures  11  and  13.  All  of  this  information  can  be 

represented  in  a  single  plot  of  .1^ { Ve }  and  a  single  plot  of  JQ{ }  through 

proper  normal ization  of  the  terms.  The  proper  normal ization  is 
2  2 

t  J  and  x  J  as  functions  of  eV  / kT .  Graphs  of  these  two  normalized 
e  o  e 

functions  are  shown  in  Figure  14. 

The  derivation  for  these  normalized  functions  follows  from  the 
combination  of  Equations  102  and  103  to  give 

Je  =  (x2/D)2  Cc2{eVe/kT)  (119) 

The  combination  of  Equations  106  and  119  then  gives 

Jo  =  (x2/D)2  f'c2{eVe/kT}  exp  {-eVe/kT}  (120) 

The  normalized  functions  are  then  from  Equation  115. 

t2J6  =  Cc2  {eVe/kT]  (121) 

and 

t2Jq  =  Tc2  (e  Ve/kT)  exp  {-eVe/kTl  (122) 

The  function  f,  is  conveniently  calculated  from  Equations  29  through  31. 

c  2 

Figure  14  shows  the  asymptotic  value  of  6.5224  for  t  J  and  the 

2  e 

asymptotic  slope  of  4  for  t  J  as  a  function  of  eV/kT  near  V  =  0. 

The  limiting  value  of  t  is  just  the  square  of  the  limiting  value  of 

Sc  from  Equation  31,  and  the  limiting  slope  is  apparent  from  the  square 
of  Equation  29. 

Equation  112  gives  JQ  and  consequently  the  cathode  work  function 
(from  Equation  8)  based  on  the  measurement  of  current,  cathode  temperature, 
and  diode  spacing.  It  does  not  require  a  knowledge  of  the  anode  work 
function  or  the  value  of  the  applied  voltage.  Equation  112  in  principle 
allows  the  determination  of  JQ  for  any  value  of  Jg.  However,  for  large  JQ 
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the  value  of  /J~~  D/x^  has  nearly  the  limiting  value  of  -2.5539.  In  this 
case  nc  is  such  a  steeply  rising  function  that  the  method  is  not  useful. 

In  practice  the  problem  with  high  values  of  JQ  can  be  overcome  by  lowering 
the  temperature,  since  JQ  is  an  exponential  function  of  the  temperature. 

6.  COMPARISON  OF  J{ V}  WITH  CHILD'S  LAW 

For  many  practical  cathode  tests  the  temperature  is  around  1300°K, 

the  cathode-anode  spacing  is  of  the  order  of  1  mm  and  the  current  density 

2 

JQ  is  of  the  order  of  1  A/cm  requiring  100V  or  more  to  obtain  JQ  from 
the  cathode.  Then  in  Equation  91  the  term  kT/e  is  of  the  order  of  .1  volt, 
and  the  voltage  offset  due  to  the  potential  minimum,  ( kT/e ) 1 n{ J/ Jq}  ,  is 
about  0.2  volts  for  J  =  JQ/ 1 0  and  decreases  as  J  approaches  J  .  This 
term,  then,  has  very  little  effect  on  an  applied  voltage  in  the  100  volt 
range.  An  examination  of  other  terms  in  Equation  91  shows  that  D/J-  /x^ 
is  in  the  range  of  a  few  hundred.  Since  £c  has  a  maximum  magnitude  of 
2.5  and  decrease  for  J  approaching  J  ,  this  term  has  only  a  small  effect 
on  the  argument  of  na  in  Equation  91.  Then  Equation  91  is  accurately 
approximated  by 

Vapp  '  (kT/e)  VDAF  /x2}  (123) 

The  function  n  (U  for  large  f,  can  be  approximated  using  Equation  27. 

a 

For  large  r.,  nfl  is  also  large  as  seen  for  example  in  Figure  4.  For  large 
ng  Equation  76  shows  that  the  function  P { u >  is  accurately  approximated  by 

P{u}  5  1  -  exp  {-u2}  (124) 

Then  for  large  £,  Equation  27  is  accurately  approximated  by 

dn  /df,  =  (2/ir1/4)n  1/4  (125) 

a  a 
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Figure  15  gives  a  graphical  comparison  of  the  approximate  expression 
from  Equation  125  compared  to  the  exact  value  of  the  derivative  calculated 
from  Equation  27.  This  graph  shows  that  the  magnitude  of  the  error  in 
using  Equation  125  decreases  slowly  as  f  increases.  However  the  error 
expressed  as  a  percentage  of  the  value  of  dn  /dt  becomes  small  very 

a 

rapidly  for  c  greater  than  50.  In  generating  Figure  15,  Equations  32 
through  34  were  used  to  calculate  n  U).  Tnen  Equation  27  and  125  were 

a 

used  to  calculate  the  derivative  and  the  approximate  derivative  as 
functions  of  f; .  Integration  of  Equation  125  shows  that 

na(U  =  (3/4)(3M1/3  C4/3  (126) 

Here  the  integration  constant  is  zero  since  n  (0)  is  zero.  Equation  126 

a 

also  follows  from  Equation  34,  though  less  intuitively.  The  number 
1.2552  in  Equation  34  is  just  [(3/4) (3/tt) I/3]3/4  rounded  to  5  digits. 

The  combination  of  Equations  123  and  126  gives 

Vapp  =  [(9/4)(02/eo)(m/2e)1/2j'2/3  J2/3  (127) 

or,  which  is  equivalent 

J  =  (4/9)eQ(2e/m)1/2V3/2  /D2  (128) 

Equation  128  is  known  as  Child's  law  (Reference  7).  Review  of  the  above 
procedure  then  shows  that  the  term  2/n  /n  in  Equation  27  is  responsible 

d 

for  the  basic  three  halves  behavior  of  J{V}. 
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SECTION  III 
SCHOTTKY  THEORY 

The  space-charge-theory  developed  to  this  point  in  Section  II  suggests 

that  there  will  be  a  zero-field  current  density  Jo  for  some  voltage  VQ  and 

that  there  will  be  no  increase  in  J  above  J  for  V  >  V  .  However, 

0  0 

Schottky  (Reference  4)  suggests  that  a  linear  applied  accelerating  field 

across  a  diode  will  interact  with  the  intrinsic  field  near  the  cathode 

surface  to  lower  the  cathode  work  function.  Consequently,  the  current 

density  will  increase  above  J  for  accelerating  fields  at  the  cathode 

o 

surface.  A  review  of  Schottky's  derivation  is  found  in  Hutchison 
(Reference  6,  pp  305  to  309).  Although  this  section  neglects  space-charge 
effects,  as  did  Schottky,  space-charge  effects  will  be  brought  into  this 
theory  and  accounted  for  in  Section  IV. 

1.  DERIVATION  OF  THE  SCHOTTKY  EQUATION 

As  discussed  by  Ramo,  Whinnery,  and  Van  Duzier  (Reference  22),  the 
force  exerted  on  an  electron  of  charge-e  located  at  position  x  away  from  a 
conducting  plane  is  the  same  as  the  force  exerted  on  the  electron  by  a 
charge  +e  located  at  position  -x.  Then  the  electric  field  experienced  by 
the  electron  is 

F  =  e/4ir<o(2x)2  (129) 

The  work  done  on  an  electron  in  bringing  it  from  a  great  distance  away 
from  the  conducting  plane  to  position  x  is  then 

x  p 

W  =  -/  [ (+e) ( -e)/T 6ttf  (x ' )  ]dx 1  =  -e(e/167rc0x)  (130) 

CO  0  w 

If  this  work  is  defined  in  terms  of  a  motive  potential  <f>  so  that 

m 

W  =  (-e)  (131) 

then 

=  e/16TiEox  (132) 
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This  potential  is  shown  in  Figure  16  where  the  Fermi  level  is  displaced 
above  the  vacuum  level  by  an  amount  ■{■,  the  work  function.  Due  to  the 
way  which  <j>  is  defined,  the  reference  zero  potential  is  the  vacuum 
level  with  no  applied  voltage.  In  principle.  Equation  132  predicts  that 
the  motive  potential  is  unbounded  at  the  cathode  surface  so  that  the 
work  function  would  have  to  be  infinite.  However,  the  idea  of  an  image 
charge  due  to  a  perfectly  conducting  surface  breaks  down  at  distances  of 
the  order  of  interatomic  spacings.  As  discussed  by  Nottingham, 

(Reference  13),  this  changes  the  shape  of  the  motive  potential  very 
near  the  surface  and  leads  to  a  finite  work  function.  Beyond  a  few 
angstroms  from  the  surface,  the  fact  that  the  motive  potential  is 
accurately  described  by  Equation  132  has  been  proven  for  many  types  of 
cathodes  by  the  agreement  of  experimental  data  with  the  Schottky  equation 
derived  below  (Reference  5,  p  109). 


If  a  potential  V  is  applied  to  a  plane  anode  at  a  distance  D  from  a 
plane  cathode,  the  applied  potential  depresses  the  vacuum  potential  Vx 
linearly  with  distance  from  the  cathode  surface  according  to 

Vx  =  Fx  (133) 


where 


F  =  V/D 


(134) 


is  a  constant  electric  field.  The  total  field  V{ x) ,  which  is  the  sum  of 
the  fields  from  Equations  132  and  133,  is  expressed 


V { x >  =  e/16*r  cQx  +  Fx  (135) 

and  is  plotted  in  Figure  17.  The  value  of  x  for  which  V  is  minimum  is 
found  from  the  zero  of  the  derivative  to  be 

*„,n  *  (,36> 

The  value  of  V  at  x„.  is 
mi  n 

<>37) 
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As  seen  in  Figure  17,  this  represents  a  decrease  in  the  work  function 
from  the  value  with  no  applied  external  voltage.  The  reduced  work  function 
is  then 

$ r  =  <P  -  ^eF  /  /4n cQ  (138) 

Using  this  modified  work  function  in  the  Richardson  equation.  Equation  8, 
gives 

J  =  AT2  exp  {-U-/e'F/4ue  )e/kT}  (139) 

or 

J  =  JQ  exp  {(e/kT)/eF/4neo  }  (140) 

where  F  =  V/D.  Calculating  this  effect  for  5000  volts  applied  across  1 

millimeter  gives  x-„  =  84  A  and  V  .  =  .085  volts  or  J  =  2.3  J  at  900°C. 

3  min  min  oQ 

Calculations  for  1000  volts  across  1  millimeter  give  xmin  =  190  A, 

Vnn'n  =  ,038  volts  and  0  =  1.4  JQ  at  900°C. 

2.  EDGE  EFFECTS  AND  SCHOTTKY  THEORY 

One  goal  of  cathode  testing  is  to  determine  basic  physical  properties 
of  a  cathode  from  specific  electrical  measurements.  This  allows  cathode 
performance  to  be  predicted  independently  of  the  test  vehicle.  For  this 
purpose  Langmuir  space  charge  theory  (References  7  and  8)  and  Schottky 
theory  (Reference  4)  are  available  to  interpret  electrical  test  data  from 
the  cathode.  However,  these  theories  are  well  developed  only  for  certain 
simple  geometries  such  as  parallel  planes  and  concentric  cylinders  and 
spheres.  In  practice  the  test  vehicle  only  approximates  one  of  these 
ideal  geometries.  When  the  theory  is  used  tc  interpret  the  measured  data, 
errors  are  introduced.  The  calculated  "basic  properties"  of  the  cathode 
are  then  dependent  to  a  certain  extent  on  the  geometry  of  the  test  vehicle. 

A  close-spaced  diode  which  is  often  used  for  testing  cathode  emission 
is  shown  in  Figure  18.  Langmuir  and  Schottky  theories  for  infinite  planes 
are  used  to  interpret  the  measured  data.  The  most  obvious  error  comes 
from  electric  field  variations  across  the  emitting  surface  due  to  electric 
field  enhancement  by  the  edge  of  the  sleeve.  Additional  error  comes  from 
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the  finite  extent  of  the  emitting  region  which  generally  is  smaller  than 
the  total  cathode  surface  (see  Figure  18).  The  presence  of  a  boundary 
to  the  emitting  layer  is  contrary  to  the  conditions  assumed  in  the 
derivation  of  Langmuir  space  charge  theory  for  parallel  planes.  The 
space  charge  is  dissipated  more  easily  at  the  boundary  resulting  in 
greater  current  density  than  is  predicted  for  the  parallel  plane  case. 

This  problem  is  discussed  by  Pierce  (Reference  23)  who  shows  that  uniform 

emission  can  be  obtained  if  a  focus  electrode  is  added  at  cathode 

potential  which  makes  a  67.5  degree  angle  with  the  normal  to  the 

cathode  surface  (see  also  Reference  22).  With  or  without  the  focus 

electrode  the  current  varies  with  the  three  halves  power  of  voltage  under 
space  charge  limited  conditions,  as  discussed  by  Langmuir  (Reference  8, 
p  251)  and  Spangenberg  (Reference  24).  Without  the  focus  electrode,  the 
saturation  is  not  abrupt  due  to  the  nonuniform  current  density.  Also, 
the  slope  of  the  three  halves  power  curve  is  not  simply  related  to  the 
square  of  the  cathode-anode  spacing. 

Focus  electrode  geometries  other  than  the  67.5  degree  sleeve  can 
improve  the  uniformity  of  the  space  charge  limited  emission  current. 
Computer  programs  for  traveling  wave  tube  gun  design  (Reference  25)  are 
available  which  can  evaluate  more  complex  geometries.  Although  these 
designs  lead  to  a  more  accurate  representation  of  the  emission  by 
Langmuir  space-charge  theory  for  parallel  planes,  they  will  not  improve 
the  representation  by  Schottky  theory  under  temperature  limited 
condi tions . 

Schottky  theory  for  parallel  planes  becomes  more  applicable  when  an 
electrode  at  cathode  potential  (guard  ring)  extends  the  plane  of  the 
emitting  surface.  The  guard  ring  gives  a  more  uniform  electric  field 
once  the  applied  voltage  is  high  enough  to  decrease  the  space  charge 
effects.  Traveling  wave  tube  gun  design  programs  can  also  calculate  the 
electric  field  for  the  test  diode  in  the  absence  of  space  charge. 

Comparing  the  Schottky  theory  for  the  true  field  distribution  with  the 
Schottky  theory  for  parallel  planes  shows  the  errors  in  the  parallel  plane 
assumption. 
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If  the  errors  in  the  parallel  plane  assumption  are  not  large,  a 
simpler  analytical  approach  gives  an  accurate  approximation  to  the  true 
field  distribution.  This  approach  begins  with  the  calculation  of  the 
electric  field  about  a  plane  edge  above  an  infinite  plane  as  shown  in 
Figure  19.  The  Schwartz  transformation  relates  the  distance  from  the 
cathode  edge  to  the  electric  potential  and  flux  through  a  procedure 
outlined  by  Ramo,  Whinnery  and  Van  Duzier  (Reference  22). 

The  Schwartz  transformation  allows  one  to  find  a  mapping  from  the 
complex  z'  plane  into  the  complex  z  plane  which  preserves  both  the 
properties  of  the  static  electric  field  and  the  boundary  conditions.  If 
the  z'  plane  is  chosen  as  shown  in  Figure  20(a),  the  solution  to  the 
static  electric  field  problem  can  be  found  in  a  straightforward  manner. 
Let  the  boundary  condition  be  V  equals  VQ  to  the  left  of  x  =  0  and  V  =  0 
to  the  right  of  x  =  0.  Then  electric  potential  v  and  flux  function  u 
are  given  by 

u  +  iv  =  ( VQ/ tt  )  1  n  { z 1 1  (141) 

where 

u  =  flux  function  =-/Edy+/Edx  (142) 

x  y 

i  =  v'TT 

F  =  electric  field  intensity 

A  quick  check  of  Equation  141  shows  that 

v  =  (VQ/„)  arctan  {y’/x'l  (143) 


satisfies  both  Lapace's  equation  and  the  boundary  condition  specified 
above  as  required  for  the  electric  potential.  The  inverse  of  Equation  141 


where 


z'  =  exp  {  ttw/  Vq) 
w  =  u  +  iv 


(144) 

(145) 


will  be  required  below  in  the  derivation. 
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For  this  particular  problem  the  z  plane  would  be  that  of  Figure  19 
and  the  required  transformation  is  (Reference  22,  p  189) 

dz/dz'  =  K  (z'-x'1/“l/T'"1^(z'-x,2)(ct2/lT'1)  (146) 

The  angles  and  u2  and  the  positions  x^  and  x£  are  shown  in  figure  20(a) 
and  (b).  Including  the  appropriate  values  for  xj ,  x2,  and  in 
Equation  146  gives 

dz/dz'  =  (K/z ’)  /z^T  (147) 

Integration  of  Equation  147  gives 

z  =  2K[/z7- T  -  arctan  {/PCT}]  +  C  (148) 

where  C  is  the  integration  constant.  The  constant  C  can  be  evaluated  by 
the  requirement  that  the  point  z *  -  1  transform  into  the  point  z  =  ih. 

This  requirement  imposed  on  Equation  148  gives 

C  -  ih  (149) 

where  h  is  the  diode  spacing.  The  constant  K  is  evaluated  from  Eqjation 
147  through  an  integration  from  point  B  to  point  C  shown  in  Figure  20  in 
both  the  z  and  z'  planes.  It  must  be  that 

/BC  dz  =  /BC  (dz/dz1)  dz’  (150) 

If  the  integration  is  about  a  small  semicircle  in  the  z'  plane  at  B  -  C, 
then  dz/dz1  becomes  arbitrarily  close  to  i K/z ’  and  dz1  can  be  expressed 
as  iz'do1.  Then  Equation  150  can  be  expressed 

CSS  dz  =  (ik/zoiz'de-  (151) 

Evaluation  of  Equation  151  then  gives 

K  =  i  h/TT  (152) 
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Finally,  a  combination  of  Equations  152,  149,  148,  and  144  gives 

x  +  iy  =  (2ih/7i)[n/2  -  arctan{/exp{7iw/V  }’  -1}  +  /exp{  ttw/ VQ)  -1] 

(153) 

Since  the  arctangent  is  related  to  the  logarithm  by  (Reference  19,  p  80) 
arctan  x  =  (-i/2)  In  {(1  +  i x ) / ( 1  -  i x ) >  (154) 

the  real  part  of  Equation  153  can  be  expressed  along  the  cathode 
surface  as 

(x/h)  =  ( 1  / 7i )  1  n{  ( 1  +  /)  -  exp{ ttu/ VQ> ) /  ( 1  -  ST~^  exp{7ru/Vo) )  }  (155) 

-  (2/tt)  ST  -  exp{7Tu/\y 

since  w  =  u  +  io  along  the  cathode  surface.  Note  that  u  ranges  from 
0  to  -  «>  along  the  cathode  surface.  From  Equation  142  the  perpendicular 
component  of  the  electric  field  Fy  is  just  du/dx.  Taking  this 
derivative  in  Equation  155  gives 

(h/VQ)  Fy  =  d(u/VQ)/d(x/h)  =  -1/vT^Txyr^T}  (156) 

The  behavior  of  Equation  155  can  be  checked  for  u  having  large 
negative  values.  In  this  case  Equation  155  can  be  approximated  as 

x/h  =  (1/ti)  In  {4  exp{-TTu/VQ}}  -  2/tt  (157) 

which  says  that 

u  =  -(Vo/h)x  (158) 

for  large  u  and  x.  This  is  expected  since  u  is  just  the  line  integral 
of  the  electric  field. 
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The  expression  of  Equation  155  is  in  an  inconvenient  form  for 
further  use  since  distance  is  given  as  a  function  of  voltage.  Attempts 
to  get  an  analytical  inverse  have  been  unsuccessful.  However,  a 
numerical  inverse  is  strai ghtforward  by  means  of  the  Newton  iteration 
(Reference  20)  also  used  in  Subsection  5(d)  of  Section  II.  The  Newton 
iteration  is  implemented  for  this  problem  by  the  function  f{u}  for  a 
specified  x  such  that 

fful  =  (1/tt)  ln{  ( 1  +  /I  -  exp{Ttu/VolT/0  -  A  -  exp{uu/V0')  J 1 


—  2 / tt )  v7l  -  expUu/V  '}'  -  x/h 

Then  the  desired  u { x }  is  found  by  solving  for 

flu)  =  0 

by  constructing  the  sequence 

U1  =  V/h  +  1  ^  */h  > 

or 

U1  =  Vox/h  if  x/h  <  -2 

and 

un  +  1  =  un  '  f{un}  /f''un}  n  =  1*2»3 


(159) 


(160) 


(161) 


(162) 


The  choice  of  starting  value  u-j  comes  from  an  inspection  of  Equation  157. 
The  function  f { up }  =  df/du  comes  by  direct  calculation  on  Equation  159  or 
by  inspection  of  Equation  156  and  is  expressed  by 

f  { un>  =  -/I  -  exp{Tti T/Vq)  (163) 


The  series  is  truncated  when  hf{u}/x  is  sufficiently  close  to  zero  and 
this  un  is  taken  as  ulxl.  Then  Fy  {xl  is  calculated  from  Equation  156. 

If  space  charge  is  neglected  and  the  surface  electric  field  is  known, 
the  emission  current  density  can  be  calculated  from  the  Schottky 
Equation  140.  A  convenient  procedure  is  to  divide  the  cathode  surface 
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into  M  rings  of  equal  area, 
ring  into  equal  areas  is 

r-  N 

The  distance  from  the  cathode  edge  to  this  mean  radius  is 

x^Nl  =  S  +  (d/2)  (1  -  /N/M  -  .5/M)  (165) 

where 

S  =  sleeve  width 

d  =  diameter  of  the  emitting  area 

A  detailed  sketch  of  this  diode  geometry  is  shown  in  Figure  18.  The 

electric  field  F  at  a  distance  xIN;  from  the  cathode  edge  is  calculated 
y 

by  means  of  the  Newton  iteration  with  Equations  155  and  156.  This  is,  of 
course,  an  approximation  since  the  cylindrical  nature  of  the  cathode  is 
ignored  in  the  field  solution.  The  current  density  at  each  x{Nl  is  then 
calculated  from  Equation  140  and  the  average  is  taken  over  the  M  rings. 

The  process  is  repeated  for  different  values  of  applied  voltage  generating 
an  I  -  V  curve.  The  Newton  iteration  and  Equations  155  and  156  have  to  be 
used  only  once  to  solve  the  field  problem  since  the  electric  field  is 
found  in  terms  of  the  applied  voltage. 

The  solution  for  the  current  density  and  electric  field  for  a 
particular  example  are  shown  in  Figure  21.  The  parameters  chosen  are 
cathode  diameter  d  =  3  mm,  sleeve  thickness  S  =  0.175  mm,  cathode-anode 
spacing  h  -  1.5  mm,  cathode  temperature  T  =  900°C,  and  cathode-anode 
voltage  Vq  =  1800  V.  The  voltage  enhancement  is  394  at  the  edge  of  the 
emission  layer  while  the  current  enhancement  is  only  about  8 %.  This  is 
surprising  at  first  glance,  since  the  current  varies  as  an  exponential 
of  the  square  root  of  the  voltage.  The  coefficient  in  the  exponent  is 
small,  however,  so  that  the  dominant  function  is  the  square  root  rather 
than  the  exponential  for  these  experimental  conditions. 


Then  the  mean  radius  which  divides  the  Nth 
=  (d/2)*'N/M  -  .5/M  (1 
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The  average  current  as  a  function  of  voltage  in  terms  of  Schottky 
coordinates  is  shown  in  Figure  22.  The  effect  of  the  edge  is  to 
increase  slightly  the  average  current  density  at  a  given  applied  voltage. 
Although  the  effect  appears  small  in  Figure  22,  the  increase  in  slope 
of  the  1 n { J / J Q }  versus  V  curve  gives  an  apparent  decrease  in  temperature 
from  900°C  to  834°C.  The  temperature  is  calculated  from  the  Schottky 
plot  by  means  of  Equation  140  in  the  form 

T  =  /e/4*. ~  (e/k)  ST  /I n{ J/Js >  (166) 

In  this  case  the  Schottky  curve  including  edge  effects  is  very  nearly  a 
straight  line  passing  through  J  =  Jq  at  V  =  0.  The  curve  including 
edge  effects  is  straight  since  the  argument  of  the  exponential  in 
Equation  140  is  small  for  this  example.  This  gives  support  to  the  common 
practice  of  using  the  Schottky  line  intercept  to  find  the  zero  field 
current  even  when  the  slope  of  the  Schottky  line  doesn't  agree  with  the 
measured  temperature. 
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SECTION  IV 

EXTENDING  SPACE-CHARGE  THEORY  INTO  THE  ACCELERATING  FIELD  RANGE 

Figure  I  shows  schematically  the  assumption  made  in  the  derivation 
of  Equation  12.  For  this  derivation,  the  potential  due  to  the  surface 
forces  which  the  electrons  must  overcome  occurs  as  a  step  function  right 
at  the  cathode  surface  and  has  a  value  of  ;  electron  volts  below  the 
Fermi  level.  Actually,  the  potential  is  not  a  true  step  function  but 
acts  over  a  short  distance  outside  the  surface  due  to  the  image  force 
on  the  electron.  Becker  (Reference  5)  points  out  that  neglecting  the 
finite  extent  of  the  image  force  has  little  effect  on  the  space-charge 
distribution  as  long  as  there  is  a  space-charge  minimum  volts  below  <t> . 
However,  if  there  is  no  space  charge  minimum,  the  space  charge  potential 
interacts  with  the  image  force  to  lower  the  value  of  4». 


1 .  BACKGROUND 

This  interaction  of  the  image  force  with  the  space  charge  occurs  in 
much  the  same  manner  as  discussed  in  the  derivation  of  the  Schottky  effect 
in  Subsection  1  of  Section  III.  The  derivation  of  the  Schottky  effect  is 
based  on  the  superposition  of  the  image  potential  of  the  electron  and  a 
linear  applied  potential.  The  following  discussion  shows  how  the  initial 
slope  of  the  space-charge  potential  serves  as  the  linear  applied  potential 
in  the  Schottky  derivation. 

Due  to  its  own  image,  an  electron  experiences  a  surface  potential 

field 


V.  =  e/4'i  (4x) 
image  o 


067) 


where  x  is  the  distance  from  the  cathode  surface,  and  V.  is  measured 

i  md  9“ 

relative  to  the  potential  a  long  distance  from  the  cathode  (vacuum  level). 
In  addition  an  electron  experiences  a  space-charge  potential  field  given  by 
Equation  27  in  normalized  form.  Figure  23(a)  shows  the  image  potential 
as  a  function  of  distance  from  the  cathode  surface  while  Figure  23(b) 
shows  the  space-charge  potential  as  a  function  of  distance  from  an 
example  cathode.  The  example  cathode  of  curve  ©  has  a  current  density 
Jo  of  5  A/cm2  at  900nC  for  dV/dx  =  0  at  the  cathode  surface.  To  illustrate 
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the  behavior  of  the  space  charge  potential  for  accelerating  field 

conditions.  Figure  23(b)  also  shows  the  potential-distance  curve 

according  to  Equation  27  for  a  cathode  with  an  accelerating  field 
-4 

dV/dx  =  6.44  x  10  volts/meter  at  the  cathode  surface.  Although  this 

2 

value  of  field  strength  (corresponding  to  the  arbitrarily  chosen  5.5  A/cm  ) 

is  chosen  at  this  point  for  demonstration  purposes  only,  the  following 

discussion  shows  how  this  increased  field  strength  gives  an  increased 

2 

current  density  of  5.5  A/cm  at  900°C  for  the  same  example  cathode. 

Note  that  the  x  scale  in  Figure  23(a)  is  expanded  by  a  factor  of  100  over 
the  x  scale  in  Figure  23(b). 

Figure  23  shows  that  over  the  range  where  the  image  potential  is 

-  8 

effective  (from  the  cathode  surface  out  to  about  5  x  10  meters)  the 
space-charge  potential  can  be  approximated  very  accurately  as  an 
extension  of  the  initial  slope  of  the  space  charge  curve 


V 

space 

charge 


{xl 


dV/dx 


(x)=  F  x 

cathode 

surface 


(168) 


Based  on  this  assumption,  namely,  that  the  space-charge  potential  is 

linear  near  the  cathode  surface,  the  derivation  for  the  Schottky  effect 

applies  if  this  F^  instead  of  V/D  is  used.  The  result  is  a  potential 

minimum  at  a  distance  x  .  as  in  Equation  136 

min  M 

x  .  =  (e/F  )1/2/(16*r  )1/2  (169) 

min  '  x'  '  o  ' 


from  the  cathode  surface,  and  the  magnitude  of  the  potential  minimum 
is  reduced  from  the  work  function  4>  by  an  amount  A<j>  where  from 
Equation  137 


A0  =  (eFx)1/2/(4nc0)1/2 


(170) 


For  these  conditions  replacement  of  *j>  with  <|>-A<j>  in  Equation  8  gives 


J  =  AT2exp  {  -eU-(eFx/4TTc0)1/2)/kT} 


(171) 
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Equation  171  with  =  0  gives  the  zero-field  current  density 

Jq  =  AT2  exp  { -e<>/kT}  (172) 

which  is  the  Ri chardson-Dushman  equation.  With  this  expression  for  JQ, 
Equation  171  becomes 

J  =  JQ  exp{e(eFx/4TH_o)1/2/kT}  (173) 


2.  GENERATION  OF  J'VI  IN  THE  ACCELERATING  FIELD  RANGE 

The  previous  discussion  now  provides  a  framework  for  calculating 
the  voltage  needed  to  give  a  specified  current.  Let  Vg  be  the  voltage 
of  the  anode  relative  to  the  cathode.  Then,  referring  to  Figure  8, 


V 

app 


Vixl 


-  Vs  + 


9,  - 


anode 


(174) 


where  x  =  D  -  xc  at  the  anode,  D  is  the  cathode-anode  spacing,  is 

the  distance  from  the  cathode  to  the  potential  minimum,  and  <j>  and  6  are 

the  work  functions  of  the  anode  and  cathode  respectively.  For  simplicity, 

the  remainder  of  this  section  takes  =  p  and  refers  to  the  applied 

a  c 

voltage  as  V.  If  J  J  ,  Equations  27  and  28  give  the  voltage  according 
to  the  procedure  outlined  by  Langmuir  with  n  =  0.  Equation  91  summarizes 
this  procedure. 


If  J  >  J  ,  Equations  27  and  20  give  the  voltage 

V  -  (kT/e)  na  (O/T  /x2}  (175) 

where  n  is  calculated  by  a  mixed  Runge-Kutta-Taylor  method  using 

a 

nQ  =  (e/kT)(x2/v(J)  Fx  (176) 

In  Equation  27.  The  expression  of  fi0  in  Equation  176  is  just  an 
expression  of  the  electric  field  in  normalized  units  using  Equations  87,  21, 
and  20.  From  Equations  173  and  176 

nQ  =  (41r,o3/2/e)(kT/e)7/4(e/21rm)1/4  ln2fJ/JQ}//J  (177) 

or 

h0  =  6.4976  x  10'2  T(7/4)  ln2i J/Jq}  //J 
with  T  in  degrees  Kelvin  and  J  in  amperes  per  square  meter. 


(178) 
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Figure  24  shows  an  example  current- vol tage  plot  for  J  ranging  from 

2 

small  values  up  through  J>JQ.  For  this  example  JQ  is  5  A/cm  ,  T  is  900°C 

and  the  cathode-anode  spacing  D  is  1  mm.  Figure  24  also  shows  the 
Child's  law  line.  Equation  128,  and  the  Schottky  line.  Equation  140  for 
comparison  with  this  extended  theory.  For  this  example  the  extended  theory 
matches  Child's  law  closely  up  to  J  =  1.15  JQ  but  the  approach  to  the 
Schottky  line  is  only  very  gradual. 

3.  BEHAVIOR  OF  J{ V 1  NEAR  J  EQUAL  TO  AND  GREATER  THAN  J 

o 

At  the  point  where  the  space-charge  range  meets  the  accelerating- 
field  range,  the  space-charge  minimum  occurs  right  at  the  cathode.  The 
derivative  of  Equation  175  with  respect  to  J  shows  the  behavior  of  the 
current-vol tage  curve  from  the  accelerating-field  side  at  the  transition 
from  the  space-charge  range  to  the  accelerating-field  range.  Since  n 
is  also  a  function  of  J,  the  derivative  is 


dV/dJ  =  (kT/e)[(D/2x,  /J  )d  n  /df, 

l.  a 


f  ' 

a 


+  (d  n  / d J ) d  n  /d  h 

O  d  U 


where 


f  ‘  =  /J  0/ x0 

a  l 


]  (179) 


(180) 


Inspection  of  Equation  177  shows  that  dnQ/dJ  =  0  at  J  =  JQ.  In  addition 
an  expression  of  n  as  an  integral  evaluated  at  r'  along  with  the 
definition  of  the  derivative  shows  that  do  /do  is  zero  for  o  =0 

3  0  0 

(i.e.,  at  J  =  J  ).  This  expression  for  do,/d5  is 

0  3  0 

.  i 

do  /d5  =  f  [n  ' 2  +  0  -  P{*^1)  explol  -  0  -  24A)]1/2dC  (181) 

a  o  aoQ  0  0 

The  derivative  evaluated  at  oQ  =  0  is  then 


do  /  dr'i 

a  o 


o=0 


q [  +  (1  -  exp  In)  -  (1  -  2v/o/iT)]1/,2dC 

'o  o 


-  /  [(1  -  Pf/^1)  explol-  (1  -  2/^)]1/2dc]/o 


(182) 
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For  the  left  hand  integral  in  Equation  182  choose  same  £  =  A  so  that 
riff}  satisfies 

\  =  (1/2)[(1  -  Pi/Ti)  exp  ini  -  (1  -  2/^h)]  (183) 

Then  a  portion  of  the  left  hand  integral  can  be  bounded  by 

/0[i1q  +  (1  -  PWn")  )exp { n  1  -  ( 1  -  2/n/V )  ^  df,  <  JQ [Hq ] ^  ^ ^  d£  =  hQA  (184) 

Observe  that  the  expression  /l  +  x  is  bounded  by  1  +  x/2  for  positive  x 
since 

/T+-*  <  .'1  +  x  +  (x2/4)  =  1  +  x/2  (185) 


Use  of  this  expression  in  the  remaining  portion  of  the  left  hand  integral 
in  Equation  182  gives 

r  1  f  1 

f  [-ig+d-Pi  vT7)  )exp{r, } -  ( l -2v/777V)  ]1/2dr.  <  /  [( 1  -Pi  /Tx)  )expf»,}-(  1  -2* nh)  ]1/2d> 

A  0 

.  \ 

f 

+  1/2/  (r^/[(l-P{,n])exp|.,}-(l-2v^A)]1/2)d£  (186) 

o 

where  the  integrals  on  the  right  hand  side  have  been  extended  to  £  =  0. 
This  preserves  the  inequality  since  the  integrals  are  positive.  Then 
from  Equations  182,  184,  and  186  the  derivative  is  bounded  by 


dn 


a 


lim  1  r-  a  + 

n  ~  In  A  + 

n  0  n  o 
o  o 


r.  2  r,' 

-j-  /  [(l-P{/^l)exp{rl}-(l-2v/v7)]"1/2d£] 
0 


V° 


(187) 


A  comparison  of  Figure  4  and  Equation  183  shows  that  A  must  go  to  zero  as 
n  goes  to  zero.  Then,  since  the  dependence  of  n  on  h  has  to  be 

0  a  0 

greater  than  or  equal  to  zero.  Equation  187  shows  that 


dn, /dn 
a  o 


=  0 


V° 


(188) 
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A  combination  of  observations  on  Equations  177,  179,  and  188  shows 
that  the  derivative  of  voltage  with  current  at  the  zero  field  current 

Jo  1s 


dV/dJ 


J=Jo 


(kT/e)(D/2x2/j)dna/df 


r  =  A 57  D/x0 

o  2 


(189) 


which  has  a  finite,  non-zero,  value.  Since  the  slope  at  J  =  JQ  is  zero 
from  the  space-charge  side  and  it  has  a  finite,  non-zero  value  from  the 
accelerating  field  side,  this  theory  predicts  a  discontinuity  in  the 
slope  of  J  as  a  function  of  V  at  J  =  J  . 

Figure  25  shows  an  example  where  this  discontinuity  in  the  slope  is 
apparent.  The  dashed  lines  intersecting  at  J  =  JQ  in  Figure  25  show 
the  limiting  slopes  as  J  approaches  Jq  from  the  space-charge  range  and 
from  the  accelerating-field  range.  Van  der  Ziel  (Reference  9)  treated 
this  example  in  the  space-charge  range,  and  his  Figure  1  shows  J 
approaching  Jy  with  a  zero  slope  as  does  this  Figure  25.  Equation  189 
gives  the  limiting  slope  as  J  approaches  from  the  accelerating-field 
range.  Figure  25  shows  that  the  limiting  slope  from  the  accelerating 
field  range  is  almost  the  same  as  the  slope  of  the  Child's  law  line  at 
J  =  J  ,  although  Child's  law  does  not  approximate  J  as  a  function  of  V 
for  this  example.  The  following  discussion  shows  that  the  agreement 
between  the  limiting  slope  and  the  Child's  law  slope  is  not  a  coincidence. 

In  many  practical  examples  such  as  the  one  treated  in  Figure  24, 
the  discontinuity  in  the  slope  of  J { V }  causes  only  a  small  perturbation 
in  the  curve.  In  this  case  the  extended  theory  provides  a  nearly  smooth 
transition  from  the  space-charge  range  into  the  accelerating  field  range 
and  then  produces  a  gradual  approach  to  the  Schottky  line.  In  this  case 
the  term  dtc/dn  in  Equation  92  is  small  except  for  a  very  small  region 
near  =  0.  Then  Equation  92  reduces  to  Equation  189  and  remains 
accurate  for  J  somewhat  greater  than  J  .  For  many  experimental  conditions 
including  those  of  Figure  24,  Equation  27  which  gives  dn./dc  ,  is 
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accurately  represented  by  only  one  term  -  the  one  involving  2/n TA  • 

a 

Consequently,  Equation  189  reduces  to 

dV/dJ  =  (kT/e)(D/2x2/J)  2  (naM1/4  (190) 

and  n,  in  Equation  190  with  their  equivalent  expressions 

a 

and  87  removes  the  normalizing  factors  and  gives 
dV/dJ  =  (m/2e)1/4t0  '1/2DV1/4/J1/2  (191) 

Rearrangement  and  integration  of  Equation  191  gives 

J  =  4/9  ,  /2e/m  V3/2/D2  (192) 

which  is  the  same  as  Equation  128,  Child's  law.  This  exercise  shows 

that  Child’s  law  is  a  good  approximation  for  the  current  as  a  function  of 

voltage  not  only  in  the  space-charge  range,  but  also  a  good  approximation 

for  a  certain  interval  into  the  accelerating-field  range.  Figure  24 

hows,  for  that  example,  that  J  follows  Child's  law  closely  up  to  about 

J  =  1.15  J  . 

o 

4.  FIELDS  AND  CHARGE  IN  THE  INTERELECTRODE  SPACE 

With  the  theory  developed  in  Subsection  1  of  Section  II  and  Subsection  1 
of  Section  IV,  it  is  now  possible  to  plot  the  voltage  V,  the  electric 
field  F  and  the  number-density  of  electrons  N  as  a  function  of  distance 
from  the  cathode  for  various  conditions  of  cathode  loading  J/ JQ .  At  a 
specified  temperature  T  and  current  density  JQ,  a  family  of  these  plots 
can  be  made  with  d/JQ  as  a  parameter.  These  curves  show  the  changing 
conditions  in  the  interelectrode  space  as  the  current  drawn  from  the 
cathode  changes,  and  they  help  to  develop  a  firmer  concept  of  the  fields 
and  charge  in  the  interelectrode  space. 

For  space-charge-limited  condition'  Equations  28  gives  the  distance  as 
a  function  of  voltage  to  the  cathode  side  of  the  potential  minimum.  After 
the  addition  of  the  normalizing  constants  from  Equations  20,  21,  and  87 
this  relationship  can  be  expressed  as 

x  =  (x2//J)  f,c{eV/kTl  (193) 


Replacement  of  x2 
from  Equations  20 
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From  Equations  84  and  85,  V  ranges  from  0  to  V$  in  Equation  193  causing 

x  to  vary  from  0  to  xm  with 
J  m 

xm  =  f'Ctln{VJ}}  (194) 

and 

V$  =  (kT/e)ln{Jo/Ji  (195) 

The  terms  x  and  V  are  measured  relative  to  the  potential  minimum; 
although  it  is  more  desirable  to  have  x'  and  V'  measured  relative  to  the 
cathode  surface  where 

x'  =  x  -  xm  (196) 

and 

V'  =  V  -  Vs  (197) 

Note  that  x  is  negative  and  V  is  positive  while  x'  is  positive  and  V' 
is  negative  in  the  cathode  space.  Then  Equation  193  can  be  expressed  as 

x'  =  (x2//j)  cc<e(V'  +  V s ) / kT }  -  xm  (198) 

The  electric  field  F  is  then  -dV/dx  where  from  Equation  28 

dV/dx  =  - (kT/e) (/J/ x9) ( [1+P{ /rO ]  exp{n  1 

1  C  C  (199) 

-0  *  3(nc/,)1/2])1/2 

with 

nc  =  eV/kT  (200) 

The  number  density  of  electrons  from  Equation  19  is 

N  =  (E()/e)  d2V/dx2  (201) 

and  from  Equation  23  this  is 

N  =  (tQkT/2e2)(J/x22)  exp{nc)[l  +  P{/^1]  (202) 
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To  the  anode  side  of  the  potential  minimum  Equation  27  gives  the 
voltage  as  a  function  of  distance.  After  addition  of  the  normalizing 
constants  this  relationship  is  expressed  as 

V  =  (kT/e)n  (x«tf  /x2> 
or 

V'  =  (kT/e )n  f(x‘  +  x  )/T  /x9l  -V 
a  m  2  s 

The  electric  field  is  then  -dV/dx  where 

dV/dx  =  (kT/e)  (v'J  /xj([i  -  PirTT  f]exp{nJ-[l  -  2/rA)1/2 

c.  da  a 

The  number  density  of  electrons  from  Equations  23  and  201  is  then 

N  =  (,0kT/2e2)(J/x22)  expfna)[l  -  Pf/^  3]  (206) 

For  accelerating  field  conditions  Equation  27  gives  the  voltage  as 
a  function  of  distance  with  the  added  condition  on  hQ  from  Equation  177. 
The  expression  for  the  voltage  is  then 

V  =  (kT/e)ii  { x/J  /x9!  (207) 

a  C 

wi  th 

nQ  =  (4 -7e)>o3/2  (kT/e)7/4(e/2T,m)1/4  ln2{J/JQl//j  (208) 

Then  from  Equation  27  the  electric  field  F  is  given  by  -dV/dx  where 

1  /2 

dV/dx  =  (kT/e) (/J  /x„)(n2  +  [l  -  Pl/O]  exp{nj-[l  -  2(n  aMV2])  (209) 

C  U  a  a  a 

In  addition,  from  Equations  23  and  201,  the  number  density  of  electrons 
is  given  by 

N  =  (!okT/2e2)(J/x22)  explain  -  P(/^  1] 
with  the  added  condition  on  h  in  Equation  208. 


(203) 

(204) 

(205) 


(210) 
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Figures  26  through  28  show  the  fields  and  charge  in  the  interelectrode 

2 

space  for  an  example  cathode  at  900°C  with  a  current  density  of  5A/cm  . 
Figure  26  shows  the  electric  potential  as  a  function  of  distance  away 
from  the  cathode  surface.  The  curve  at  J/JQ  =  .6  shows  clearly  the 
potential  minimum.  The  curve  at  J/JQ  =  1  shows  the  zero  slope  of  V { x } 
at  x  =  0.  The  curves  for  J/JQ>1  show  the  linear  behavior  of  V { x }  near 
x  =  0.  The  function  J { V }  can,  in  principle,  be  constructed  from 
Figure  26  by  drawing  a  line  of  constant  x  at  the  distance  x  =  D  represent¬ 
ing  the  diode  spacing.  The  intersection  of  this  line  with  the  J { V }  lines 
then  gives  the  voltage  associated  with  each  value  of  J  for  the  specified 
diode  spacing.  This  line  of  constant  x  must  fall  to  the  right  of  the 
potential  minimum  to  generate  a  valid  J { V }  curve,  however,  since  the 
retarding  field  range  requires  a  non-zero  h0  in  the  solution  for  Tin  } . 

In  contrast  to  the  accelerating  field  field  range,  this  A  for  the 
retarding  field  range  depends  on  the  spacing  D  and  cannot  be  incorporated 
into  the  series  of  curves  in  Figure  26  which  are  independent  of  D 
(see  Reference  9). 

Figure  27  shows  the  derivative  dV/dx  which  is  the  negative  of  the 
electric  field  F.  The  curve  at  J/ J  =  .6  shows  -F  to  be  negative  near 
the  cathode.  The  electric  field  goes  to  zero  at  the  potential  minimum, 
and  then  becomes  positive  in  the  anode  space.  The  curve  for  J/JQ  =  1 
shows  F  -  0  at  x  -  0  as  required  for  the  zero  field  current  density. 

As  J  increases  above  JQ  the  curve  for  -F  rises  rapidly  near  the  cathode 
and  then  levels  off  at  distances  away  from  the  cathode.  This  tendency 
becomes  more  pronounced  for  large  J/JQ  until  the  electric  field  is  nearly 
constant  in  the  interelectrode  space.  This  constant  electric  field  is, 
of  course,  expected  as  space-charge  effects  become  small. 

Figure  28  shows  the  number  density  N  of  electrons  as  a  function  of 
distance  from  the  cathode.  For  small  J/JQ  the  number  density  of  electrons 
is  large  near  the  cathode  due  to  electrons  returning  from  the  potential 
minimum.  However  the  number  density  falls  off  rapidly  away  from  the 
cathode  toward  the  potential  minimum.  For  J/Jq  closer  to  1,  the  number 
density  of  electrons  is  smaller  near  the  cathode  but  it  does  not  fall  off 
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so  rapidly  going  away  from  the  cathode.  From  Equations  195,  200,  and  202 
the  electron  number  density  at  the  cathode  is 

Nix  =  01  =  («  okT/2e2)(Jo/x2Z)[l  +  Pi lnt JQ/J> } ]  (211) 


At  the  zero  field  current  density  J  ,  Equation  211  shows  that 


Nix  =  01|  =  (i:okT/2e2)(Jo/x22) 

J=J 


(212) 


As  J/Jq  becomes  very  small.  Equation  211  shows  that  the  number  density 
of  electrons  becomes  twice  that  given  in  Equation  212, 


Nix  =  01 


-  (, 


0kT/e2)(J0/x22)  = 


2N( x  =  01 


J  =  0 


(213) 


J=J. 


since  the  error  function  Plu!  tends  to  1  for  large  argument.  Equation  213 
is  reasonable  since  at  J  -  J  none  of  the  electrons  are  returning  while 
at  J  =  0  all  the  electrons  are  returning  to  give  twice  the  density  at  the 
cathode  surface.  As  J  becomes  larger  than  the  electron  density  at  the 
cathode  is  again  larger  than  it  is  at  J  =  and  it  falls  away  much  more 

rapidly  at  distances  away  from  the  cathode.  From  Equation  210  the  electron 
density  at  the  cathode  surface  is 

N ' x  =  01  =  (- QkT/2e2)(J/x22)  (214) 

which  increases  linearly  with  increasing  J. 

Figure  28  suggests  that  the  electron  density  decreases  away  from  the 
cathode  with  a  finite  slope  for  J  less  than  or  equal  to  JQ.  However  for 
J>J0  the  electron  density  has  a  decreasing  slope  with  unbounded  magnitude 
at  the  cathode  surface.  This  can  be  shown  rigorously  from  the  equations 
involved.  From  Equations  28  and  202 

dN/dx  =  -(>0kT/2e2)(J3/2/x23)[exp{nc}  (1  +  Pl/^  })  (215) 

+  1/ */-nc’][(l  +  P  1 1  ) exp {  n c } - ( 1  +  2/n^n ) ^ 
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For  large  n  this  tends  to  the  limiting  value 

dN/dx  =  -(/2  cokT/e2)(Jo3/2/x23)  (216) 

which  gives  the  slope  at  the  cathode  surface  for  J  much  less  than  JQ. 

Here  use  is  made  of  Eguations  195  and  200  and  it  is  noted  that 

l//nc  and  /n^/expln  }  both  become  small  for  large  n  .  It  is  also  assumed 

for  this  illustration  that  the  anode  is  sufficiently  far  away  to  allow  a 

potential  minimum  in  the  cathode  space.  For  nc  near  zero,  corresponding 

to  J  near  J  ,  Equations  45,  54,  and  215  show  that  dN/dx  has  the  limiting 

value 

dN/dx  =  - ( cQ kT/2 /n  e2)  ( JQ3/2/x23)  (217) 


This  can  also  be  shown  directly  from  Equation  215  by  the  use  of  L. 
Hospital's  rule  to  show  that 


1  im 


[1  +  P{/n7l]exp{n  }-[l  +  2/n/Ti  ] 

0  (I/2/77")  (  - -c  -  -C - 2 - 


1/2 

)  =  (1/2/7) 


(218) 


In  all  these  cases,  then,  from  very  small  J  up  to  J  =  JQ,  the  slope  of 
the  electron  density  dN/dx  is  finite  at  the  cathode. 

For  J>Jq  the  derivative  dN/dx  is  found  from  Equations  64  and  210 

dN/dx  =  (.  okT/2e2)(J3/2/x23)[(l  -  P{/^  })  exp  lna)  (219) 

-1  //tth"  ][h  ^  +  (1  -  P  t/nTl)  exp  In) 

a  u  a  a 

-1  (  1  -  2/ri  / 71  )  ]^2 

a 

The  term  l/>nj  increases  without  bound  as  na  approaches  zero  showing 
that  the  slope  of  the  electron  density  is  unbounded  at  the  cathode. 

This  is  true  for  any  fiQ>0  corresponding  to  J>dQ- 
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5.  EFFECTS  OF  DIODE  SPACING  ON  J ■ V ( 

The  zero-field  current  density  J  is  often  estimated  from  experi- 

o 

mental  data  through  a  determination  of  the  "space-charge  breakaway  point." 
This  is  the  point  where  the  JiVl  curve  departs  by  a  specified  amount 
from  Child's  law.  One  method  of  specifying  this  "breakaway  point"  uses 
the  projections  of  the  slopes  of  the  JiVs  curve  from  the  space-charge 
region  and  from  the  accelerating-field  region  plotted  on  two-thirds 
power  paper.  The  intersection  of  the  two  lines  representing  these 
projected  slopes  is  taken  as  an  estimate  of  J  .  The  dashed  lines  in 
Figure  29  give  an  example  of  this  method.  Affleck  (Reference  17) 
presents  experimental  data  showing  that  the  diode  spacing  affects  the 
value  of  the  "breakaway  point"  determined  in  this  manner.  The  following 
discussion  uses  the  extended  theory  developed  in  this  report  to  show  how 
the  "breakaway-point"  technique  leads  to  erroneous  values  of  Jq  which  are 
dependent  on  the  diode  spacing. 

The  extended  theory  of  Subsection  2  of  Section  IV  is  used  to  calculate 

the  J ( V •  curve  for  a  cathode  similar  to  the  one  tested  by  Affleck. 

Figure  29  shows  the  current  density  as  a  function  of  applied  voltage  for 

two  values  of  diode  spacing,  0.3  mm  and  0.9  mm.  For  the  example  the 

2 

temperature  is  900 "C  and  the  zero-field  current  density  is  5  A/cm  .  The 

dashed  lines  in  Figure  29  show  that  the  "space-charge  breakaway  point" 

gives  values  of  "J  "  higher  than  the  true  value  of  J  in  each  case.  The 
J  o  J  2  0 

projected  slopes  intersect  at  "J  "  -  7.6  A/crn  for  the  0.3  mm  spacing  and 
2  0 

at  "Jo"  =  6.3  A/cm  for  the  0.9  mm  spacing.  This  property  of  the 
"breakaway  point"  giving  higher  values  of  "Jo"  for  the  close-spaced  diodes 
than  for  those  with  wider  spacings  agrees  with  Affleck’s  conclusions 
based  on  data  from  a  diode  with  variable  electrode  spacing. 

The  above  example  shows  that  the  "breakaway  point"  gives  only  a 
very  crude  estimate  of  the  value  for  J  .  In  the  0.3  mm  case  the  estimate 
is  50%  too  high,  while  in  the  0.9  mm  case  it  is  25%  too  high.  One  reason 
for  the  erroneously  high  values  of  JQ  is  that  J { V >  follows  the  Child's  law 
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line  to  higher  values  of  J  above  JQ  in  the  close-spaced  diode  than  for 
the  one  with  a  wider  spacing.  A  second  reason  is  that  the  comparison 
for  the  cases  with  different  diode  spacing  is  made  over  the  same  range  of 
diode  voltage.  Since  the  electric  field  is  higher  in  the  close-spaced 
diode  the  current  density  rises  to  a  considerably  higher  values  in  this 
case.  This  causes  the  line  projecting  the  slope  of  J1V)  from  the 
accelerating-field  range  to  lie  in  the  region  of  higher  values  of  J  for 
the  close-spaced  diode.  In  conclusion,  the  use  of  intercepting  slopes 
on  two-thirds  power  paper  to  determine  the  "space-charge  breakaway  point" 
leads  to  erroneous  values  of  JQ  which  are  dependent  on  the  diode  spacing. 
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SECTION  V 
CONCLUSIONS 

A  method  for  treating  space-charge  effects  in  the  accelerating-field 
range  is  given.  The  method  results  in  the  current  as  a  function  of 
voltage  for  close-spaced,  paral lei -plane  diodes.  The  space-charge  theory 
developed  by  Langmuir  is  applied  to  the  accelerating-field  range  along 
with  the  image-potential  theory  suggested  by  Schottky.  The  key  to 
combining  these  two  theories  is  to  note  that  the  space-charge  potential 
is  essentially  a  linear  function  of  distance  over  the  range  of  action  of 
the  image  potential.  The  combined  theory  produces  a  current-vol tage 
curve  which  begins  at  the  zero-field  current  and  voltage  predicted  by 
space-charge  theory  and  then  gradually  approaches  the  Schottky  line. 

For  those  conditions  where  Child's  law  approximates  the  current-vol tage 
curve  in  the  space-charge  range.  Child's  law  also  approximates  the 
current-vol tage  curve  generated  by  the  combined  theory  for  an  interval 
of  current  values  extending  above  the  zero-field  current. 

The  solution  procedure  outlined  by  Langmuir  is  summarized  in  a  single 
equation.  The  derivative  of  this  equation  shows  more  explicitly  the 
behavior  of  the  current-vol tage  curve  at  the  transition  points  from  the 
exponential  range  to  the  space-charge  range  and  from  the  space-charge 
range  to  the  accelerating-field  range.  The  current  and  the  derivative  of 
the  current  are  both  continuous  functions  of  voltage  at  the  transition 
from  the  exponential  range  into  the  space-charge  range.  At  the  transition 
from  the  space-charge  range  into  the  accelerating-field  range  the  current 
is  a  continuous  function  of  voltage,  but  the  derivative  of  the  current  is 
discontinuous.  This  discontinuity  arises  because  Langmuir  space-charge 
theory  always  gives  a  zero  slope  for  the  current-vol tage  curve  at  the 
zero  field  current.  However,  the  combined  theory  gives  a  finite,  non-zero 
slope  for  the  current- vol tage  curve  at  this  transition  point. 
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VACUUM 


LIMITING 
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6  REM  Data  is  dEta/dXi  at  X  =  0  ,  delta  X, print  steps, Xmox 


15 

READ  D 0  ,  H  ,  L  ,  B 

20 

DSPLAY  "DELTA  Xi 

=  "  ,H, "dY^dX (at  X  =  0) 

22 

D  S  P  L  A  Y 

30 

LET  X  =  Y  =  Q  =  0 

4  0 

LET  E  =  1  ''  L 

90 

DSPLAY  "Value  of 

X  i  "  ,  "  V  a  I u  e  of  Eta" 

1  O0 

G  0  S  U B  3  2  0 

120 

GOTO  160 

130 

LET  U  =  Y 

1  3  1 

G  0  S  U  E  10  0  0 

1  ■->  /. 

LET  F  1  =  P  5  - ( 1-2  + 

U  1  s  1  .77245) +  D  0  ■#  D  0 

136 

IF  F  1  <  0  LET  F 1  = 

0 

137 

LET  F1=SQR (F 1 ) 

1  38 

LET  U  =  Y  +  H / 3  *  F 1 

139 

G  0  S  U  B  10  0O 

140 

LET  Y=Y+H*(F1+H 

/  4  *  P  5  ) 

150 

LET  X  =  X  +  H 

160 

IF  (X+H-B) +1E-0 

7  >=  0  THEN  190 

170 

LET  R  =  X  +  H 

130 

GOTO  200 

I  9  0 

LET  R  =  B 

2O0 

LET  fl=INTIE*R)/ 

E 

2  1  O 

IF  (  A  +  L  -  R  ) -1E-0 

7  >=  0  THEN  230 

220 

LET  A  =  A  +  L 

2  30 

IF  (  X - A )  + 1 E - 0 7 

>=  0  THEN  260 

240 

LET  G  =  A 

2  50 

G  O  S  U  B  3  2  0 

2  60 

IF  R  =  E  THEN  280 

270 

GOTO  130 

280 

IF  ABS'R-AJ<.5* 

.000001  THEN  311 

290 

LET  G  =  R 

300 

G  0 S  U  E  3  2  O 

3  1  1 

STOP 

320 

LET  U  =  Y 

32  1 

G  0  S  U  B  10  0  0 

322  LET  F1=P5-(1-2*U1/1.77245)+D0»D0 

323  LET  F1=SQR(F1) 

324  LET  U  =  Y  + (  Q  -  X  ) / 3  *  F  1 

325  LET  X 1 =  X 

326  LET  X  =  Q 

327  GO SUE  1000 

3  2  8  LET  Yl=Y+{0-Xn*  (F1  +  (Q-X1)/4*P5) 

329  LET  X=X1 
370  DSPLAY  0  ,  Y 1 

330  RETURN 

390  DATA  0,. 01,  .01,  .2 

980  REM  THIS  SUBROUTINE  CALCULATES  ( 1 -ERF ( SQR C Y ) )) *EXP ( Y ) 
1000  LET  U  1  =  S  Q  R  (  U  ) 

1010  LET  T=l/fl+.327591*Ul) 

1020  LET  P5  =  T* «  1  .42141+T*  (-1  ,45315+T*1  .06141  )  ) 

1030  LET  P5=T* ( . 25483+T* ( - .2845+P5) ) 

1040  RETURN 
1040  RETURN 


Figure  6.  BASIC  Language  Program  Using  Mixed  Runge-Kutta-Taylor  Method 
to  Solve  for  n,U) . 

d 
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0  E L  T  R  xii 

Value  of  X  i 
0 

.01 
.  02 
.03 
.04 
.05 
.06 
.07 
.03 
.09 
.  1 
.  1  1 
.  1  2 
.  1  3 
.  1  4 
.  1  5 
.16 
.  1  7 
.  IS 
.  19 


.0f 

Value  of  Eta 
0 

.000025 
9 . 97337E-05 
2 . 23642E-04 
3 .96303E-04 
6 . 19072E-04 
8.89S64E-04 
1  .20892E-03 
1  .  57603E-03 
1  .  9  9  1  1  6  E  -  0  3 

2 . 4  5  3  7  6  E  -  0  3 
2 . 9  6  3  6  6  E  -  0  3 

3 . 5  2  G  6  2  E  -  0  3 
4  .  12445E-03 
4.77433E-03 
5 .47154E-03 
6 .2143 5 E- 0  3 
.007003 

?  .83726E-03 
8 .716  9OE-03 
9.641  69E-03 


d  V'7  d  X  1  a  t  X  =  01 


Figure  7.  Sample  Output  from  the  Proaram  of  Figure  5. 
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Figure  8.  Schematic  of  the  Potential  Distribution  Between  the  Anode 
and  Cathode  under  Space-Charge  Limited  Conditions. 
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JQ  =  10  itiA/cm^ 
T  =  700"C 
D  =  0.2  nr.i 


d(J/JQ)/dV 


d2(J/Jo)/dV2 


APPLIED  VOLTAGE  V 


L VJL I s  J 


Figure  9.  The  Current  Density  and  its  First  and  Second  Derivatives  as  a 
Function  of  the  Applied  Voltage  for  an  Example  Cathode, 
for  this  example  :  is  assumed  equal  to  *  .  The  dashed  lines 
cover  the  retarding  field  range  while  thecconti nuous  curves 
cover  the  spare  charge  range.  The  value  of  the  voltage  at 
.1  J  is  0.9994  only  by  coincidence  in  this  particular  example 


[A/m4] 
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x  =  4.58 
T  =  900°C 
D  =  1 .0  mm 


x  =  1.05 
T  =  700°C 
D  =  0.2  mm 


.3105 

LIMITING 

VALUE 


[A/m4-] 


5.864 
LIMITING 
VALUE 


Figure  11.  The  zero  field  current  JQ  as  a  Function  of  Jg  the  Current  at 
the  Transition  from  Retarding  Field  to  Space-Charge-Limited 
Operation  for  D  =  lmm,  T  =  900  C  and  for  [)  =  0.2  mm,  T  =  700°C. 


2  D  S  P  L  ft Y  "Program  of  3  Nov  30" 

4  DS  P  L  ft  Y  "Prints  J  <  as  a  function  of  Jo" 

5  D  S  P  L  ft  V 


10 

LET 

T  = 

10  73 

20 

LET 

D  = 

.  0  0  0  1 

30 

D  S  P  L  ft  V 

"  T  E  H  P  =  (  '  K  )  "  ,  T  ,  “  D  I  S  T  =  (  M  )  "  ,  D 

45 

LET 

ft  2 

=  CT  1.751*1. 08  8  6E-06 

49 

LET 

N  = 

.  1 

50 

D  S  P  L  ft  V 

"  D  X  2  =  "  ,  D  /  ft  2 

5  1 

D  S  P  L  ft  Y 

52 

D  S  P  L  ft  Y 

"  J  o  "  ,  "  J  <  "  ,  "  V  <  "  ,  "  X  I  c  " 

5  3 

FOR 

0  = 

1  TO  4  STEP  1 

54 

LET 

N  = 

N  *  1  O 

5  5 

FOR 

I  5 

=  N  TO  1 0  *  N  STEP  N 

60 

LET 

I  = 

6 . 5  *  i:  ft  2  -  D  i  *  (  ft  2  /"  B  ) 

65 

I  F 

I  > 

=  15  LET  I =  .5*15 

70 

LET 

E  1 

=  -  L  0  G  t  I  I  5  I 

75 

I  F 

E  1  < 

O  LET  E 1 =  O 

9  0 

LET 

V  1 

=El*T*8.62E-05 

1O0  LET  E  =  E  1 

110  G o S U  £  14  0  0 

12  0  LET  H  =  Z  +  D  " ft  2  *  $ Q R (  I  1 

2  3  0  LET  T  1  =  1  i  1  +  .327591  *  S  Q  R  (El  1  ) 

29  0  LET  P  5  =  1  .42  141+ (-1  .45315+1  ,06141*T1)*T1 
291  LET  P  5  =  (  .25  433+  i  -  .284497+P5+T1  )  *  T  1  )  *  T  1 

3  0  0  LET  F  2  =  2  ■*  E  X  P  (El  )  -  P  5  -  1  -  S  Q  R  (El  )  *  1  .  12338 

3  10  IP  F  2  <=  0  LET  F2=1E-12 

3  2  0  LET  F  2  = 1 / 3 0 R f F  2 ) 

322  LET  H  1  =  F 2 -  I  +  Bs 2  '  R 2 / S QR (  I  ) 

3  2  4  LET  I  1  =  I  - H  -  H 1 

3  2  6  IF  ft  E  S  (  t  I  -  I  1  i  r- 1)  <  .  0  0  0  1  GOTO  332 

323  LET  1=11 
3  3  0  GOTO  7  0 

3  3  2  D  S  P  L  ft V  I  5  ,  I  ,  V  1  ,  Z 

350  NEXT  15 

351  N E X T  0 
1336  STOP 

1390  REM  THIS  SUBROUTINE  CRLCULRTES  ZC  RS  fl  FUNCTION  OF  ETRC 
14  0  0  IF  E  >  =  .2  THEN  1450 

1410  LET  RsSGR(E) 

1420  LET  Z=-2*R+.376126*E+2.51956E-02*E*R+6.9163E-04*E*E 

1421  LET  Z  =  Z -  9 . 9  4  O  9  E - 0  4  *  E  *  E  *  R 
1440  RETURN 

14  5  0  IF  E  >  =  3  THEN  1490 

14  6  0  LET  Z  =  - (  -  6 .15614+E-65. 691 2«E*E-59. 489 3*Et3+4. 10335 *Et4) 
1461  LET  Z=Z/ 262909-14. 6202*E-49.5951*Et2-16.0639*Et3+Ef4) 
1480  RETURN 

1490  LET  Z=-2.55389+1.41421*EXP(-E/2)-.0123*EXP(-E) 

1491  LET  ZsZ+1/3^1 .41421*(SQR(E/3.14159)+l)*EXP(-3*E/2) 

,rio  RC  7 

i6<ro 

Figure  12.  BASIC  Program  to  Solve  Equation  113  and  Generate  Jg  and  Vg 
as  Functions  of  JQ. 
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,e  1  .  I  he  /pro  field  Current  .1  and  the  Transition  Current  Jc  Plotted 

.1  .  turn  lions  of  the  Transition  Voltage  V  for  0  0.?mm,  T  -  700  C. 

I 


ZERO-FIELD  CURRENT  [ A/ni' 


r.Ve/kT  [NORMALIZED  UNITS] 

Figure  14.  Graphs  of  Jp  and  Jp  as  Functions  of  V  .  These  Universal  Curves 
Relate  the  Measured  Current  at  the  Exponential  Point  to  the 
Zero  field  Current  of  the  Cathode. 
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external  field  is  10'  volts/meter. 
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Figure  18.  Schematic  of  a  Close-Spaced  Diode  Tester  Showing  the 

Cross  Section  of  a  Cylindrical  Cathode  and  Anode.  The 
diameter  of  the  emitting  layer  is  the  same  as  the  sleeve 
inner  diameter. 
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Figure  19.  A  Plane  Edge  Above  a  Conducting  Plane.  The  flux  function  u 
is  set  to  zero  at  X  =  0,  y  =  h. 
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Figure  21.  Current  Density  and  Electric  Field  at  the  Cathode  Surface  as 
a  Function  of  Distance  from  the  Cathode  Edge  -  All  in 
Normalized  Units.  Here  S  =  sleeve  thickness,  h  =  cathode  to 
anode  spacing,  VQ  =  applied  voltage,  T  =  temperature  of  the 

cathode,  d  =  diameter  of  the  emitting  area,  X  =  distance  from 
cathode  edge,  and  J  =  zero-field  current  density. 
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Figure  22.  Theoretical  Schottky  Plots.  The  curve  "with  edge  effects" 
takes  into  account  the  enhanced  field  at  the  cathode  edge. 
The  curve  "without  edge  effects"  refers  to  infinite  parallel 
planes. 
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Figure  25.  Current- Vol tage  Plot  Showing  a  Break  in  the  Slope  of  J  as  a 
Function  of  V  at  J  =  J  .  The  dashed  lines  show  the  limiting 

values  of  the  slope  at  J  =  J  .  The  Child's  law  line  and  the 
Schottky  line  are  shown  for  comparison  with  the  combined  theory 
developed  in  this  paper. 
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